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ABSTRACT 



In this thesis we study the geometry of the fixed point set E of a 
smooth mapping $ : M — M on an n-dimensional Riemannian manifold 
(M, g) by computing the asymptotic expansion as t O"*" of the trace of the 
deformed heat kernel U(t; x, $(a;)) of the Laplace operator on M. We assume 
that the fixed point set is a union of connected components Sj, i = 1, . . . , s, 
each of which is a smooth compact submanifold of M with dimension rrii, 
^ rrii ^ n — 1. 

There exists the asymptotic expansion as t — 0"*" 

/s oo » 

where Cif'iy) are scalar invariants on Ej depending on covariant derivatives 
of the curvature of the metric g and symmetrized covariant derivatives of the 
differential (i$ of the mapping $, evaluated on Ej. 

Due to the localization principle, it is possible to compute the expan- 
sion for each component separately. We develop a generalized Laplace method 
for computing the coefficients a!'^ (y) in this expansion and compute the coeffi- 
cients Co, ai, and 02 explicitly in the following cases: 



1. M is a flat two-dimensional manifold, S is a zero- or one- dimensional 
component of the fixed point set. 

2. M is an n-dimensional curved manifold, E is a zero- dimensional compo- 
nent of the fixed point set. 

We also develop algorithms for computing further coefficients and 
expressing them in a symmetrization-free form. 
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CHAPTER 1 



Introduction 

1.1 Differential Geometry 

Here we provide a brief review of differential geometry concepts that 
we will use further. For more detailed treatment of the subject see, for exam- 
ple, [iniini. 

Definition. A topological Hausdorff space M with a countable base is called an 
n- dimensional manifold (without boundary), n G N, if, for any point p G M, 
there exists an open neighborhood U C M of the point p, open subset V C M", 
and homeomorphism (p : U ^ V. The pair {U, (f) is called a coordinate chart 
or local coordinates. 

We always assume that all manifolds we deal with are smooth. That 
is, given an n-dimensional manifold M and two coordinate patches {U, ip) and 
{V,ip), such that U (IV ^ 0, the corresponding transition map 

is smooth (infinitely differentiable). 

Definition. Let M be a manifold and f : M —>■ M. be a function on it. The 
function f is continuously differentiable on M if, for any point p E M 
and any coordinate patch {U, ip), containing the point p, the function f o ip~^ : 
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(p{U) ^ M. is continuously dijferentiable. If f is continuously differentiable on 
M, we will write f E C^{M). 

In the same manner one can give definitions of /c-times differentiable 
functions C^{M) and infinitely differentiable functions C°°{M), as well as dif- 
ferentiable mappings / : M — M". We will denote the space of smooth func- 
tions with compact support by C^{M). 

Definition. Let M be an n- dimensional manifold. A vector X. at a point 
p E M is a linear differential operator of the first order X : C^{M) R. In 
local coordinates it can be written as 

n 

X(/) = ^X^ 

i=l 

where X* are components of X. The vector X is also called the tangent or 
contravariant vector. The set of all tangent vectors to the manifold M at 
the point p is called the tangent space and is denoted by TpM. The union 
of all tangent spaces to the manifold M at all points is called the tangent 
bundle and is denoted by TM. A vector field (or a section of the tangent 
bundle) v on the manifold M is a mapping v : M — TM, such that, for any 
point p G M, there holds v(p) G TpM. 

Below we follow the standard Einstein convention of summation over 
repeated indices from one to the dimension of the manifold. Where it is neces- 
sary, we will specify limits of the summation explicitly. In local coordinates a 
vector field v may be written as 

v(p) = v\p)di, 
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with di being an i-ih. vector of the coordinate basis in the tangent space TpM, 
defined by 

where are local coordinates. We will say that the vector field v is smooth 
and will write v G C°°{TM), if all its components v^{x) are smooth functions 
in all coordinate patches. 

Remark. One can show that, although we use local coordinates to introduce 
various geometrical objects, they have invariant geometrical meaning which 
does not depend on the choice of a coordinate patch. All our results will be 
formulated in the invariant form. 

Definition. Let M be a manifold. The cotangent space T*M at a point 
p & M is the dual space to TpM, i.e. the space of all linear functionals on TpM . 
Its elements are called cotangent vectors, covariant vectors or covectors. 
The union of all cotangent spaces at all points is called the cotangent bundle 
and is denoted by T*M. A covector field (also called a 1-form or a section 
of the cotangent bundle) v on the manifold M is a mapping v : M — * T*M , 
such that, for any point p G M, there holds v(p) G TpM. 

The basis in a cotangent space, dual to the coordinate basis di, is 
denoted by dx"^ and is defined by the condition dx^{dj) = 5*, where 5* is the 
Kronecker symbol. As well as for vector fields, we define the set of all smooth 
covector fields C°^(T*M). 

Definition. Let ^ : M ^ M be a smooth mapping on a manifold M . Let 
p E M and p' = Let be local coordinates in the point p, x^' be local 
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coordinates in the point p' , and be the coordinate functions of the mapping 
$. The differential of the mapping $ at the point p is a linear mapping 



: TpM Tp,M 



with components 

The tensor product of two functionals / : X — » M and : F — > M 
is a functional f®g: X^Y^M. defined by (/ ® g){x, y) = f{x)g{y) for any 
X G X and y ^Y. 

Definition. A tensor T of the type [s, q) at a point p of a manifold M is 
a linear functional 

T : TpM X ... X TpM x T*M x . . . x T*M M . 

' V ' ' 

9 s 

In the coordinate basis it is given by 

T = T^;;if^ dx^^ ® dx^' ® <9ii ® . . . ® .9^, . 

The set of all tensors of the type [s, q) at all points of the manifold M is called 
the {s,q) -tensor bundle and is denoted by T^M. A smooth {s,q) -tensor 
field is a smooth mapping assigning to each point p E M a tensor of the type 
{s,q). The set of all smooth {s,q) -tens or fields is denoted by C°°{T^M). 

Wfien s or g in tlie above definition equals zero, we will not write it, 
that is, = TP and = T^- 
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For any m G N let Sm. be the set of all permutations of the numbers 
1, ... ,771. A transposition is a permutation, which exchanges only two numbers, 
leaving all others fixed. Every permutation can be represented as a composition 
of transpositions. Such a representation is not unique, but the number of 
transpositions in it is an invariant mod 2. The signature of a permutation 
(p G Sm, denoted by sign{ip), is 1, if it is represented by a composition of even 
number of transpositions, and —1 otherwise. For any permutation ip G Sm. and 
any tensor T G we define the tensor (pT G by 



We will say that T G T^M is a symmetric tensor, if ipT = T for 
any ip G Sm- Let VmM C T^M be the space of symmetric (0, m)-tensors. We 
define the operator of symmetrization Sym : T^M — s> VmM by 



Sym(T) = i- ^T. 
m! ^-^ 

We will also use parentheses to indicate symmetrization: 

(Sym(T))i,.„i„ = %.„i„). 

If it is necessary to exclude some indices from symmetrization, we will delimit 
them by | . . . | , that is. 



'^ih---ik\ji---jl\ik+l---im) I -^(i,p(i)---V(fe)bi 



We will say that T G TmM is an antisymmetric tensor, if ipT = 

sign{{p)T for any ip G 5*^. Let A^M C T^M be the space of antisymmetric 
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(0, m)-tensors (also called m- forms). It follows from the definition, that, for 
any n- dimensional manifold M, only spaces A"*M with ^ m ^ n are non- 
trivial, that is, they contain non-zero tensors. We define the operator of 
antisymmetrization (alternation) Alt : T^M — * A™M by 

Alt(T) = i- sign((/^)</^T. 
m! ^-^ 

We will also use square brackets to indicate alternation: 

(Alt(T)),,...,„ 

The exterior or wedge product A : A^M x A'^M A^'+'^M is defined by 

a A (3= ^^t?^' Alt(a®/3). 
p\ ql 

When it is not confusing, we will write 7^^, Vn, and A" instead of 
TgPM, VnM, and A"M respectively. 

We say that a manifold M is equipped with a Riemannian metric 
g a g is a smooth symmetric tensor field of the type (0, 2) and it defines a 
positive definite non-degenerate inner product in each tangent space TpM. We 
will call the pair (M, g) a Riemannian manifold. On a Riemannian manifold 
we can use the metric for raising and lowering tensor indices: 

tti = gija\ a' = g^'ui, 

where g^^ are components of the inverse matrix to gij, i.e. g^^gik = S-j,. The 
determinant of the matrix gij we will denote by detg or by \g\. 

For integration over a manifold M we need to define a volume form. 
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Definition. An n- dimensional manifold M is orientable, if there exists a 
constant sign n-form a G A". 

In case of an orientable Riemannian ra-dimensional manifold (M, g) 
there is a special volume form generated by the metric. It is called the Rie- 
mannian volume form and is given by 

dvol = e{di, . . . , (9„) 1^1^/^ dx^A...A rfx", 

where e{di, . . . , dn) is equal to 1 for positively oriented basis di and —1 other- 
wise. The volume form changes sign whenever a basis changes the orientation; 
such a geometrical object is called a pseudoform. We will always work with 
orientable manifolds and will use only positively oriented bases without speci- 
fying it explicitly. In this case e{di, . . . , dn) = 1 and we will omit it henceforth. 

When computing integrals over a manifold, it is usually necessary to 
work in several coordinate charts. For using formulas, written in local coordi- 
nates, a partition of unity over a manifold is used for splitting an integral into 
parts, such that each part can be covered by a single chart. 

Definition. Let (M, g) be a Riemannian manifold. The Christoffel symbols 

T^ji^ are defined by 

The covariant derivative (also called the connection, compatible with 
the metric, or Levi-Civita connection) is an operator V : C'^{T^) 
C^{Tg\,), defined by 

s q 

("TTT^yi-.-is _ fj rpii...is _|_ \ ^ -pi; ^ll ...«(_ i mij^i .. .l^ _ \ ^ -pm rpii...is 

1=1 1=1 



For a given point p & M and a vector X G TpM, the covariant derivative in 
the direction of the vector X is 

(VxT)-- =X^(VT)--,. 
We will denote covariant derivatives by indices after a semicolon, 

ijiii...is _ /^rrTi\ii...is 

jl...jq;kl...kr ~ ) jl...jqkl...kr' 

Definition. The components of the Riemann curvature tensor are defined 



by 

pi _ _ jfe I -pi -pm _ -pi -pm 



The Riemann curvature tensor has the following properties: 

Rijkl = —Rjikl = —Rijlk = Rklij , (1-1) 
R^jkl + R\lj + R\jk = , (1.2) 
R^''kl;m + R^'' lm;kR^'' mk;l = . (1.3) 

The last two equations are called Bianci identities. For the commutator of 
covariant derivatives there holds 

s q 
r\rpi-i...is \ ^ -pii rpii...ii-imii^i...is \ ^ pm rpii...is /i r\ 

1=1 1=1 

Contraction of two indices of the Riemann curvature tensor gives the Ricci 
curvature tensor 

R jik Rjk ; 

and contraction of all indices gives the scalar curvature 

oil' dA: p 

Wij — K k — ^- 
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Definition. A geodesic on a Riemannian manifold (M, g) is a critical curve 

1 

7 : [0, 1] — > M of the length functional l{^) = j dt (7(t),7(t)) . 



One can show that geodesies satisfy the equation 

V^(t)7(t) = 0, 

which in local coordinates takes the form 

d'^i'' ■ di^ 

Since we consider smooth manifolds, the theory of ordinary differential equa- 
tions guarantees that, for any point p G M and any vector v G TpM, there 
exists an e > 0, such that equation (jl.5|) has a unique solution 7v(t), t G {—e, e), 
satisfying initial conditions 7v(0) = p and 7v(0) = v. On a smooth compact 
manifold without boundary this solution exists for any t G M. 

Definition. Let (M, g) be a smooth compact Riemannian manifold without 
boundary. For any point p ^ M the exponential map exp^ : TpM M is 
defined by 

expp(v) = 7v(l), 

where 7v(^) is the solution of the equation (|1.5p with the initial conditions 
7v(0) = p and 7(,(0) = v. 

For non-compact manifolds the exponential map is well defined only 
in some neighborhood of G TpM. 
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1.2 World Function 

The local geometry of a smooth Riemannian manifold (M, g) can be 
described by the world function (see, for example, PH])- The world function is a 
real-valued non-negative function a defined in a neighborhood of the diagonal 
of M X M as follows. Let a;' G M be a fixed point. We fix a sufficiently 
small neighborhood of x' in M so that every point x in this neighborhood 
can be connected by a single geodesic to the point x', which means that the 
exponential map is injective within this neighborhood. Let r(x') be the geodesic 
radius of such a neighborhood. The injectivity radius of the manifold M, 
denoted by rinj(M), is defined as the infimum of r(x') over the whole manifold: 

rinj(M) = inf r{x). 

We assume that the injectivity radius is strictly positive, i.e. r^^j > 0. For 
smooth compact manifolds this is always the case. For non-compact manifolds, 
this will be one of our assumptions (among other assumptions specified below). 
The world function a{x, x') is defined as half the square of the length of the 
geodesic connecting the points x and x'. 

Let y = y{s) be a parametrization of the geodesic connecting the 

points X and x' with s being the natural parameter (arc length) along the 

geodesic. Let y(0) = x and ?/(r) = x'. Then y"(s) = —;—{s) is the unit 

as 

tangent vector at the point y{t) of the geodesic. The world function satisfies 
the equations [IH] 

a{x, x) = -g ^ ) ^(^' ^ ) ' 
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and 



(1.6) 



It is easy to see that 



exp^(r?/"(0)) = x' . 



(1.7) 



If two points X, x' G M cannot be connected by a single geodesic, it 



is still possible to define the distance between x and x' as the infimum of the 
length functional over all smooth curves connecting x and x', if they belong 
to the same connected component. The distance between two submanifolds 
A, B G M we define as 



1.3 The Laplacian on a Riemannian Manifold 

Now we will describe the construction of the Laplace operator on a 
Riemannian manifold {M,g). For details see [3 13 EI- 

We define the Hilbert space of square integrable functions L'^{M) to 
be the completion of the space C^{M) of smooth functions with compact 
support in the norm induced by the L^-inner product 



dist{A,B) 





M 



The metric g defines the inner product in every tangent space TpM 
and induces in a natural way the inner product in the dual space T*M which 
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we will also denote by g. It enables one to define the L^-inner product in the 
space C^{T*M) of smooth sections of the cotangent bundle T*M with compact 
support by 

{a,p) = J dvo\{x) g{a,p) 

AI 

and consider a completion of C^{T*M) in the norm induced by this inner 
product. We will denote this completion by L'^(T*M). 

The covariant derivative on the space of smooth functions is the map- 
ping V : C°°(M) ^ C°°(T*M). Since C^iM) C C°°(M) and C^iM) is dense 
in L'^{M), the operator V is densely defined. Thus, we can define the adjoint 
operator V* by 

(a,V/) = (V*a,/), 

where a G C^{T*M) and / G C^{M). 

Now we define the Laplace operator on C^(M) by 

A = -V*V. 

It is a symmetric operator, i.e. {Aj\h) = {f,Ah) for any f,h ^ C^{M). 
Moreover, it is essentially self-adjoint, which means that there is a unique self- 
adjoint extension of A. To simplify notation, we will denote this extension by 
the same symbol. In local coordinates the expression for the Laplacian is 

Af = \gr'/' d, l^g'^ \g\'^' dj) . 

Notice, that the leading symbol of the Laplacian is negative definite: 

aLiA■x,O = -9''ix)^^^,<0 



13 



for any x G M and any ^ 7^ 0. Therefore, the leading symbol of the operator 

F = -A 

is positive definite. 

Definition. We will call a differential operator of order m over an n-dimensio- 
nal manifold M without boundary elliptic, if there exist constants Ci, C2 > 0, 
such that for any x E M and any ^ G its leading symbol aiix,^) satisfies 
the inequality 

Ciier^K(x,oi^c'2ier, 

where = + ■ ■ ■ + C- 

Let H he a. Hilbert space and A : H H he a linear operator on 
it. Let A G C. If the inverse of the operator A — X exists, we will call it 
the resolvent and denote it by Rx = {A — A)~^. There are the following 
possibilities P3] : 

1. Rx exists, is bounded, and is densely defined in H. Then A is called a 
regular value of the operator A. The set of all regular values is called 
the resolvent set and is denoted by g{A). 

2. Rx does not exist. That is, the kernel of the operator A — A is non-trivial. 
Ex = Ker(y4 — A) 7^ {0}. Then A is called an eigenvalue of the operator 
A. The vector space Ex is called the eigenspace corresponding to the 
eigenvalue A, and any non-zero vector ip E Ex is called an eigenvector, 
corresponding to the eigenvalue A. If the space Ex is finite-dimensional, 
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then d{X) = dim Ex is called the multiplicity of A. The set of all eigen- 
values is called the point (or discrete) spectrum and is denoted by 

3. -Ra exists and is densely defined in H, but is unbounded. The set of all 
such A is called the continuous spectrum and is denoted by crdA). 

4. Rx exists, but is not densely defined in H. The set of all such A is called 
the residual spectrum and is denoted by ar{A). 

The spectrum cr{A) of the operator A is the union of discrete, continuous, 
and residual spectra, 

a{A) = ap{A) U (Tc{A) U ar{A) . 

Clearly, the spectrum is the complement of the resolvent set, cr{A) = C \ g{A). 
If if is a finite dimensional space, then all points of the spectrum of the operator 
A are eigenvalues, that is, c^A) = ap{A), but it may not be the case if H is 
infinite-dimensional . 

Both operators A and F are elliptic. There is the following well-known 
theorem about the spectrum of an elliptic self-adjoint differential operator with 
a positive definite leading symbol [TT] . 

Theorem 1.1. Let {M,g) be a compact Riemannian manifold without bound- 
ary and F : C°°(M) C°°{M) be an elliptic self-adjoint differential operator 
with a positive definite leading symbol. Then: 

1. The spectrum of the operator F is real, discrete, and bounded from below. 
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2. All eigenspaces of the operator F are finite- dimensional. 

3. Eigenfunctions of the operator F are smooth and form an orthonormal 
basis in L'^{M). 

More generally, a second-order elliptic differential operator acting on 
smooth sections of a vector bundle is called a Laplace type operator, if it 
has a scalar leading symbol. 

1.4 Heat Kernel and Spectral Functions 

In this section we will describe the heat kernel and other spectral 
functions on a Riemannian manifold {M,g). For details see \^ I13| IT71 IT^. 

The heat kernel U (t; x, x') of a Laplace type operator F on an n- 
dimensional Riemannian manifold without boundary M is the integral kernel 
of the heat semigroup operator 



and can be defined also as the fundamental solution of the heat equation 



where 6{x,x') is the Dirac distribution. Near the diagonal set of M x M the 
heat kernel has the form ^ |2] 



exp(-tF) : L2(M) ^ ^^(M) , 



{dt + F)U{t;x,x') = 



for t > 0, with the initial condition 



f/(0+;x,x') 



6{x, x') 



U{t-x,x') = {4TTty''/^A'/\x,x')exp 
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where A{x,x') is the Van Vleck-Morette determinant, 

A{x,x') = \g{x)\-'^'det{-V,VM^,x')) \g{xT'^' , (1.9) 
and the function Q{t] x, x') has the asymptotic expansion as t ^ 0"*" 

nit;x,x') r^Y,—^ ^k{x.x') . (1.10) 

fc=0 

(Asymptotic equivalence will be defined in Section IT31 ) 

For a compact manifold M the heat semigroup operator exp(— tF) is 
a bounded trace-class operator on the Hilbert space L'^{M), with the trace 

Q{t) = Trexp(-tF) , 

given by 



dvol(x)f/(t;x,x) = 

M ™=1 



oo 

tXr] 



e 



where each eigenvalue Am, m G N, is repeated as many times as its multiplicity. 
This trace is also called the trace of the heat kernel and is one of the spectral 
functions of the operator F. Some other spectral functions may be expressed 
in terms of the trace of the heat kernel by means of integral transforms. 

The distribution (or counting) function N{X) is defined as the 
number of eigenvalues below A G M 

iV(A) = |{Am : Am ^ A}| , 

and is given by 

e+ioo 
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for A 7^ Am,, m G N, where e is a positive constant and 6' is a Heaviside 
distribution. For A = Am we have 

iV(Am) = hm^ N{X). 

The density function g{X) is defined as the derivative of the distri- 
bution function 

^.(A) = ^iV(A) 

and is given by 

e+ioo 

^(A)= / — e*^e(t) = 5^5(A-Am). 

m=l 



Let A G C be such that Re A < Ai. Then the operator F — A is 
positive. Let s G C be such that Re s > — (this condition is necessary for 
convergence of the following integral and sum). The generalized (^-function 
C(s, A) is defined as the trace of a complex power of the operator F — X 

as,x) = TT{F-xr 

and is given by 

C(5,A): 



-A)' 

Ti 

The generalized (^-function is an analytical function of s for Re s > — and can 
be continued analytically to a meromorphic function in C In particular, it is 
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analytic at s = and, thus, it is possible to define the functional determi- 
nant of the operator F — A by 



Det(F - A) = exp (-dsC ) . 

\ s=0/ 



These spectral functions are very useful in studying the spectrum of 
the operator F and, if known exactly, they determine the spectrum. This is not 
valid for their asymptotic expansions and there are examples of operators with 
the same asymptotic series of spectral functions, but different spectra [7]. In 
principle, all these functions are equivalent to each other, but the heat kernel 
is more convenient for practical purposes, since it is a smooth function, while 
the distribution and density functions are singular. 

1.5 Laplace Method 

In this work we will develop a generalized Laplace method for com- 
puting asymptotics of integrals over manifolds. Here we present necessary 
information for both one- and multi-dimensional Laplace methods. For proofs 
of cited theorems and lemmas see, for example, |31 |S] . 

Definition. Functions f{t) and g{t) are asymptotically equivalent as t ^ 

0^ , denoted by f g, if 

t^0+ t^ 

for any N > 0. 

Definition. Let S E C°^{[a, b]). A point ^ E (a, b) is called a non- degenerate 
critical point of the function S, if S'{^) = and S"{^) ^ 0. 
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Theorem 1.2. Let (f,S & C°°{[a,b]). Let the function S attain its global 
minimum on [a, b] at an interior point C, G (a, b) . Let the point C, be the only 
point where the global minimum is attained, that is, S{x) > S{^) for any x ^ C,, 
and assume that ^ is a non- degenerate critical point. Then there exists an 
asymptotic expansion as t 0^ 

b ^ 
I dx^(x)exp {-\s{x)) ~ exp 

The coefficients Ck depend polynomially on derivatives of the function (p at the 
point C,, on derivatives of the function S of order higher than two at the point 
t and on [S"{Or'^'- 

For computing the coefficients Ck in Theorem II .21 the foUowing resuh 
for standard Gaussian integrals is useful. 

Lemma 1.1. Let a > 0. Then for any A; G Z, k ^ 0, there holds 

dxe^'^'^'x^''^' = 0, 

oo 
oo 

- „-l/2-fc 

A theorem similar to Theorem II . 21 holds in the multi-dimensional case. 

Definition. Let Q CW' be a bounded open set. Let S G C°°(r2). A point ^ G 
is called a non- degenerate critical point of the function S , if VS{^) = 
and the Hessian matrix of the function S at the point C,, that is, the matrix with 
entries {didjS){^), is either strictly positive or strictly negative definite. 
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Theorem 1.3. Let Q C MJ^ be a bounded open set. Let (p,S ^ C°°{Q). Let 
the function S attain its global minimum on Q at a point C,- Let the point ^ 
be the only point where the global minimum is attained and assume that C, is 
a non- degenerate critical point. Then there exists an asymptotic expansion as 
t^0+ 

I dx^(x)exp [--S{x)) ~ exp (-^^(o) f^Cfcr/'+'. 

The coefficients Ck depend polynomially on derivatives of the function ip at the 
point C,, on derivatives of the function S of order higher than two at the point 
C,, and on the inverse of the Hessian matrix Q = [{didj S) {^)~\ . 

Lemma 1.2. Let A be a strictly positive symmetric matrix and G = A~^. Then 
for any k & k ^ 0, there holds 

j dx exp [-x^Ax) x'^ ... x'^'^+i = , 

dx exp {-x^Ax) x'^ ... x'^^ = ^^)^J' (det Gf'^ G^'^'^ . . . ^^2.-1*2^) _ 

1.6 Spectral Asymptotics of Riemannian Submanifolds 

In this section we wiU introduce deformed spectral functions, follow- 
ing [HI [m Hn] , and formulate the problem considered in this work. 

Let (M, g) be a smooth compact Riemannian manifold without bound- 
ary, $ : M — > M be an isometric mapping on it, and S be the fixed point set 
of the mapping $, that is, 

S = {x e M : <l>(x) = x}. (Lll) 
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The mapping $ naturally defines the operator $ : C°°{M) —* C°°{M) by 

($/) (x) = /($(x)). 

The composition of the operator $ and the heat semigroup operator defines an 
operator 

U^{t) = ^ exp (tA) (1.12) 

whose kernel is 

U^{t;x,x') = U{t■,^x),x'), (1.13) 

which we will call the deformed heat kernel. 
The deformed heat trace is 

e$(t) = Tr U^{t) = J dvol(x) U{t; ^x),x) . (1.14) 

M 

It can be used to define other deformed spectral functions of the Laplace op- 
erator on the manifold M. Namely, the deformed distribution function 

£+ioo 
E—ioc 

the deformed density function 

e+ioo 
e—ioo 

the deformed (^-function 

oo 
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and the deformed functional determinant of the operator F — X 



Det$(F — A) = exp ( —dsC<s> 



s=0 



} 



where we have used the same notation as in Section IT!^ 



These deformed spectral functions were considered in [TH]. The fol- 



lowing theorem is known for the asymptotic expansion of the deformed 
heat trace. 

Theorem 1.4. Let M be a smooth compact Riemannian manifold without 
boundary. Let ^ : M ^ M be an isometry on M. Let S be the fixed point set 
of^. Then: 

1) is a disjoined union of connected submanifolds Sj, z = 1, . . . , s; 

2) ast^ 0+ 



where rrii = dimSj and % (y) are scalar invariants on Sj. 

Donnelly [HI computed the coefficients qq and ai of this asymptotic ex- 
pansion. Asymptotic expansions of this kind are used for proofs of the Atiyah- 
Singer-Lefschetz formulas for compact group actions and index theorems. 



etry we will consider a smooth mapping $ (not necessarily an isometry) with 
the fixed point set E, satisfying the following assumptions: 




s oo 



'i 



Our goal is to consider a more general problem. Instead of the isom- 
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1. The fixed point set S consists of finitely many disjoint connected compo- 
nents Sj, z = 1, . . . , s. 

2. Each connected component Sj is a smooth compact submanifold of M 
without boundary. 

3. Each connected component Sj is a "non-degenerate" fixed point set. We 
will formulate this assumption precisely for all considered cases in Sec- 
tion O (p- EEI), SectionHOfp.HI^. and Section O (p- ESI) • 

A theorem analogous to Theorem 11.41 holds for such mappings as well ^1] . All 
these assumptions are automatically satisfied for isometries. 

We will develop a technique for computation of the asymptotic ex- 
pansion coefficients and compute explicitly the coefficients oq, ai, and a2. 

We will say that two subsets A and B of the manifold M are separated 
if dist(A,i?) > 0. It follows from the above assumptions, that all components 
Sj are separated, that is, 

mindist (Sj, Sj) > 0. 

Since the asymptotic expansion of the deformed heat trace is determined 

by a neighborhood of S and all connected components Sj are isolated, it is pos- 
sible to consider each component separately. The whole asymptotic expansion 
is equal to the sum of asymptotic expansions generated by each component. 
So, without loss of generality and in order to simplify notation, we will assume 
that S consists of only one connected component. 

We will also consider non-compact manifolds M. In this case our 
additional assumption will be that the injectivity radius of M is non-zero. 



CHAPTER 2 



Symmetric Tensors 

In this chapter we will derive reduction formulas for contraction of 
symmetric tensors, that allow us to express multi-dimensional Gaussian in- 
tegrals in a symmetrization-free form. Our main result in this chapter is 
Lemma [2. 11 which gives an explicit expression of this kind for a general case. 

2.1 Algebra of Symmetric Tensors 

Let \/ : Vp X Vq ^p+q be the symmetrized tensor product of 
symmetric covariant tensors, defined by 

Ay B = Sjm{A ® B) . 

Let Q be a symmetric (2, 0)-tensor. Let r : Vm+2 ^ Kn be the 
operator of contraction of symmetric covariant tensors with Q, defined by 

= g"^A„^,,..,,„ . (2.1) 

By definition, tA = for any A from Vq or Vi. 

A mult i- index n of the order is a fc-tuple (ni, . . . , n^) of non- 
negative integers. We will use the following notation: 

|n| = ni + n2 + . . . + Uk , 
n! = rill 77,2! . . . nfc! . 
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2.2 Contraction Formulas: Theory 

Let /c G N and n be a multi-index of the order k, such that |n| = 2N. 
Let A^^^ G Vrn, i = 1, . . . ,k, be symmetric tensors. We will be interested in 
explicit formulas for the full contraction of the symmetrized product of these 
tensors with the tensor Q, that is, in scalar expressions of the form 



V . . . V A^*^)) . (2.2) 



By "explicit formulas" we mean expressions without symmetrization of indices. 
We will develop a technique for obtaining such formulas and then apply it to 
some particular cases that will be needed in the following chapters. 

By definition, we have 

r^(A«V...VA('=)) = -l^ J2 TM, (2.3) 

where 

T(^) = ( g ® . . . ® Q ® . . . ® A«) 

N 

We need to combine all similar terms in the sum ()2.H|) . that arise because of 
the symmetry of the tensors A^*) and Q, as well as the possibility to permute 
factors in the tensor product of Q's. 

With each term T{ip) we will associate a graph F = H{ip). Every 
graph r has k distinct vertices of degrees n^, corresponding to the tensors A^^\ 
and edges, corresponding to the factors of the tensor Q. An edge between 
vertices i and j means that an index of the tensor A^^^ is contracted with an 
index of the tensor A^^^ by one of the factors Q. Let 0(n) be the set of all such 
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graphs. For every permutation (/? G S2N the term T{(p) has a unique graph 
r = H{(p) associated with it. Thus, the mapping 

H:S2N^ 0(n) 
is well defined. From the definition of the mapping H it follows: 

1. Given a graph F G (S(n) it is possible to construct many terms T{ip) or, 
in other words, to find many permutations ip G S2N, such that H{lp) = F. 

2. For any graph F G ©(n) all permutations (p G i/~^(F) generate similar 
terms T((y9) in the sum (j2.3p . 

3. For any two different graphs Fi,F2 G C5(n) and any two permutations 
(pi G ff~"'^(Fi), (p2 G i/~"'^(F2) terms T((y9i) and T{if2) are not similar. 

So, the mapping if is surjective, but not injective. Therefore, we have a parti- 
tion of S2N into equivalence classes 

^27v= U H-\T). 

rGl»(n) 

All permutations 9? from the same class H~^{r) generate similar terms T[ip). 

Our goal is now to list all elements of C5(n) and to count how many 
similar terms T{ip) in the sum ()2.3p correspond to each graph. That is, for any 
graph F G C5(n) we want to find the number of elements in the equivalence 
class H^^(r). We will denote it by 

c{T) = \H-\T)\ . 
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We will identify a graph F G C5(n) with its adjacency matrix F = {'jij), 
which is a square matrix of order k. By definition, its entry '-fij is equal to the 
number of edges between the vertices i and j, that is, to the number of indices 
of the tensor A^^^ contracted with indices of the tensor A^^^ by factors of Q. It 
is easy to see, that every matrix F G (J5(n) has the following properties: 

1. For any i and j the entry 'jij is a non- negative integer from to 2N. 

2. The matrix F is symmetric, that is, for any i and j there holds = 7jj. 

3. For any i the diagonal coefficient is even. 

4. For any i the sum of all entries of the i-th row, as well as the sum of all 
entries of the i-th column, is equal to the rank of the tensor A^^^: 

k k 

Conversely, every matrix of the order k having all these properties is an element 
of the set (5(n). 

Let F G 0(n). For computing c(F) we select an arbitrary element 
ifr G H~^{r) and the corresponding term T{ip). We can generate new terms 
associated with F by using symmetry of the tensors Q and A^^\ and permuting 
factors in the tensor product of Q's. It gives us 2^A^! ni! . . . rifc! terms. However, 
some terms will differ just by names of dummy indices. Below we compute how 
many of these generated terms are the same (up to names of dummy indices). 

1. If ^ 2, then simultaneous permutation of two indices of the contracted 
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pair in the tensors Q and A*^*) does not change the term. For example: 



It reduces the number of terms by a factor of 2'^"'^. Note, that — is the 

number of self-loops at the vertex i, and - tr F is the number of self-loops 
in the whole graph F. 

2. If ^ 4, then simultaneous permutation of contracting factors of Q 
and contracted pairs of indices of the tensor A^"^^ does not change the 
term. For example: 

^ ^ ap"fa... ^ ^ yoap... 

It reduces the number of terms by a factor of 



2 

3. If ^ 2 for i < j, then simultaneous permutation of indices of the 
tensors A^^^ and A*^-'^ contracted with the tensor Q, and factors of Q, does 
not change the term. For example: 

It reduces the number of terms by a factor of 7jj!. 

Taking all this into account, the number of terms associated with F is 

2^N\ni\...nk\ 



c(F) 



n n (?)!2(^-)/^ 



Note, that these coefficients have the following property: 

J2 c{T) = {2N)\. 
re©(n) 
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So, it is natural to define weights or probabilities of graphs T by 

c(r) n! 1 



c(r) 



{2N)\ (|n|-l)!! JJ JJ 



and they have the property 

re©(n) 

Thus, we have proved the following lemma. 

Lemma 2.1. Let k,N E N. Let n be a multi-index of the order k, such that 
|n| = 2N. Let A^^^ E Vn^, i = 1, . . . , k, be symmetric {0,ni) -tensors. Then 

r^(A« V...VA(^') = Yl c(r)T(^r), 

reis(n) 

where (pr is an arbitrary permutation from the set H~^(r) and the coefficients 
c(r) are given by 

n! 1 



c(r) 



(i-i - n n 7. 



2.3 Contraction Formulas: Particular Cases 

Now we will apply Lemma 12.11 to a few particular cases. Matrices in 
the following corollaries are listed in the "decreasing row-wise lexicographical 
order" . 

Corollary 2.1. Let A,B e V3. Then 

r3(AVi?) = ^(3T(i)+2T(2)), 
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where T(^i) are terms associated with the adjacency matrices r(j) listed below: 

r(i) = (i 2)' ^(2) = (3 o)- 

Corollary 2.2. Let A e V3, B e V5. Then 

r^(AVfi) = i(3T(i)+4T(2)), 
where T(^i) are terms associated with the adjacency matrices r(j) listed below: 

r(i) = (i 4), ^(2) =(3 2)- 

Corollary 2.3. Let A e V4, B e V^. Then 

t\A V 5) = ^ (3T(i) + 24T(2) + 8T(3)) , 
where Ti^i) are terms associated with the adjacency matrices r(j) listed below: 



r(i) = (o 4), r(2)=(2 2)' ^(3) =(4 

Corollary 2.4. Le^ e ^3, C G ^4. Then 

(A V 5 V C) = ^ (3 T(i) + 12 T(2) + 12 T(3) + 8 T(4) + 2 T(5) 

+ 24 T(6) + 12 T(7) + 24 T(8) + 8 T(9) 

where T(^i) are terms associated with the adjacency matrices r(j) listed below: 




(2 1 0~ 
r(2) = I 1 2 I , 1(3) = 
VO 2 2, 



^0 3 0' 

3 

,0 4, 

^0 1 2' 

1 2 

.2 2 0, 



1(5) = 3 , 1(6) 



r(8) =1 2 , 1(9) 
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Corollary 2.5. Let A,B,C,D e ^3. Then 



r*^(A V 5 V C V D) = ^ (9 Td) + 6 r(2) + 18 T(3) + 18 r(4) + 36 T(5) 
+ 18 T(6) + 18 T(7) + 9 T(8) + 18 T(9) + 18 T(io) + 36 T(n) + 6 T(i2) 
+ 9 T(i3) + 18 T(i4) + 18 T(i5) + 6 T(i6) + 36 T(i7) + 18 T(i8) + 6 T(i9) 
+ 4 T(2o) + 36 r(2i) + 18 r(22) + 36 T(23) + 18 T(24) + 36 T^25) + 36 T^26) 

(2.4) 

+ 18 T(27) + 36 T(28) + 36 T(29) + 18 T(3o) + 36 T(3i) + 36 T(32) + 144 T(33) 

+ 36 T(34) + 18 r(35) + 36 T(36) + 36 T(37) + 6 T(38) + 4 T(39) + 18 ^^q) 
+ 36 r(4i) + 36 r(42) + 36 T(43) + 18 T(44) + 36 r(45) + 6 T(46) + 4 T(47)j , 

where T(^i) are terms associated with the adjacency matrices r(j) listed below: 
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We will also need formulas for contractions of equal tensors. 
Corollary 2.6. Let AEV3. Then 

T%A V A V A V A) = ^ (9T(i) + 12T(2) + 72T(3) + 24T(4) 

+ 144 T(5) + 4 T(2o) + 72 T(23) + 48 Tog)) , (2.5) 

where T(^i) are defined as in Corollary 



Proof. When some or all tensors in the contraction ()2.4|1 are equal, some adja- 
cency matrices generate similar terms. Namely, if tensors A^^^ and A^^^ in (|2.2p 
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are equal, and matrices r(i) and r(2) can be obtained one from another by 
exchanging i-th and j-th rows, and i-th and j-th columns, then r(i) and r(2) 
generate similar terms. By direct comparison, the matrices r(j) from Corol- 
lary 12.51 split into eight groups, such that all matrices in each group generate 
similar terms. Numbers i for groups of equivalent matrices r(i) are listed below: 

Group 1: 1, 8, 13; 

Group 2: 2, 12, 16, 19, 38, 46; 

Group 3: 3, 6, 9, 10, 14, 18, 22, 24, 27, 35, 40, 44; 

Group 4: 4, 7, 15, 30; 

Group 5: 5, 11, 17, 21, 26, 28, 31, 32, 34, 37, 41, 43; 

Group 6: 20, 39, 47; 

Group 7: 23, 25, 29, 36, 42, 45; 

Group 8: 33. 

By adding the coefficients of the terms T(j) in the formula (|2.4p separately for 
each group listed above, we obtain the coefficients in ()2.5|) . □ 



In the same way one can obtain the following corollary. 
Corollary 2.7. Let A e V3, B e V4. Then 

= (3 Td) + 24 r(2) + 12 To) + 16 T(4) + 2 T(5) + 24 T(6) + 24 r(8)) , 
where T(^i) are defined as in Corollary \2.Ji\ 



CHAPTER 3 



Spectral Geometry of Submanifolds: Flat Manifold 

In this chapter we consider the plane, M = M^. We provide new 
proofs for asymptotic expansion Theorems 13.11 and 13.21 These proofs will give 
us a way to compute the asymptotic expansion coefficients. The coefficients 
ai, 02 in Lemma (3.31 and the coefficients Ai, A2 in Lemma 13.41 are computed 
explicitly for the first time. 

The injectivity radius of the plane is infinite, rinj(M^) = 00, and the 
world function is defined globally, a : x ^ M. We will use the Euclidean 
metric in and local coordinates x". In Cartesian coordinates the world 
function is given by 

a(a;,a;') = ^5„^(x'"-x")(x'^-x^), 
and the Laplacian by 

The heat kernel U{t; x, x') has the well-known form 

U{t;x,x') = {4nt)-'exp(^-^^^^ . 

Let $ : ]R2 ^ ^ smooth mapping and let S, U<^{t), and 

X, x') be defined as in Section Oby fTTT]) . (fO^ . ^TWf . Let 5 : m 
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be a real- valued function defined by 

S(x) = cr . 

Since the world function a is defined globally, the function S is well-defined. 
Obviously, S{x) ^ for any x G and S{x) = on S. 

Since dimM^ = 2, there are only two possible cases: dimS = 1, that 
is, S is a curve; and dimS = 0, that is, S is a point. We will consider these 
two cases separately. 

3.1 One-Dimensional Submanifolds 

Let dimS = 1. Let h be a smooth unit normal vector field along the 
curve S and Vh be the corresponding normal derivative. We define functions 
Sfc : S ^ M by 

Sfc = V^;^ , A;eNU{0}. (3.1) 



Since S 



= and Vh'S' 

s 



s 

= 0, it follows that sq and si vanish. 

s 

2 . 10)2 



Theorem 3.1. Let ^ be a smooth mapping $ : M . Let S be its fixed 

point set: S = {x G : $(a;) = x}. Let S be a real-valued function S* : — > M 
defined by S{x) = a {^{x),x) . Let $ satisfy the assumptions: 

($.1) is a smooth compact connected one- dimensional submanifold of 
without boundary. 

($.2) There exist constants p > 0, R > 0, C > 0, such that for any x G M^, if 
a [x, 0) > — , then 

S{x)^C[a{x,0)]r (3.2) 
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($.3) S2 = V^'S' > at all points ofT,, where Sk are defined as in (j3.1|) . 
Then there exists the asymptotic expansion as t ^ 

oo 
fc=0 

where the coefficients Ak are locally computable in the form 

Ak= dyol{y) ak{y), 



and ak are scalars on S depending polynomially on S2 ^^"^ and Sk, k ^ 3. 

Remark. Assumptions ($.1)-(<I>.3) mean that the fixed point set S is a smooth 
closed curve without self-intersections, the function S{x) grows sufficiently fast 
at infinity and it is non-degenerate near E. 

Proof of Theorem \3.1\ Let be Cartesian coordinates in M^. We have 
e$(t) = (int)-^ J dx exp (^-^S{x)^ , 

where dx = dx^dx"^ is the standard Lebesgue measure. The proof of the theorem 
will be based on a generalized Laplace method and involves the following steps: 

1. Reducing the integration in the deformed heat trace Q<s,{t) to an integral 
over a suitable neighborhood of E. 

2. Choosing a special tangent-normal coordinate system in Mg. 

3. Scaling the normal variable by t^^"^. 

4. Expanding the integrand in the power series in t^/^. 

5. Computing integrals over the normal variable. 
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3.1.1 Tangent and Normal Coordinates 

Let the fixed point set E be described by the equations 



S : = ny), 



dx 



where y is a parameter on S. The vector field e = — is tangent to E. In a 

dy 

sufficiently small neighborhood of E we can introduce local coordinates for 
in the following way. We construct a family of geodesies x = x{y, z) in M^, such 
that at each point of E the tangent vector to the geodesic is equal to h, the 
smooth unit vector field, normal to E. We choose the parameter 2; in a natural 
way, so that it is the signed distance to E along the geodesic. That is, we have 
the differential equation 



with the initial conditions 



a;" 
dx"" 



dz 



2=0 



2=0 



= 

h%y) . 



(3.3) 

(3.4) 
(3.5) 



We restrict ourselves to a tubular neighborhood of E with |2;| < e 
for sufficiently small e, so that the change of coordinates (x^,x^) — > {y^z) is 
not degenerate. We will specify the exact condition later. The solution of the 
equations ()3.3|) - ()3.5p is given by 

x^{y,z) = r{y) + zh^{y). 



At every point in we have: 

dx°' 



dy 



e" + 2- 



dy ' 



dx"' 
dz 
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Due to the choice of we also have: 



dx"" dx^ ^ ^ dx"" dx^ 
^ dz dz ' dy dz 



dx dx 

i.e. the vectors — and — form an orthogonal system. 
ay oz 

Finally, the volume form in the coordinates (y, z) is 

dvol(?/, z) = [g{y, z)f'^ dyAdz = ^ 
= [7(2/)]'^' [l + zx{y)]dyAdz, 

where is the standard Levi-Civita symbol with ei2 = 1, g = det Qais, 7 is 
the induced metric on S and 

Now we see, that it is enough to take e less than inf \k{v)\~^ . Since S is 
compact, this infimum is positive. 



3.1.2 Reduction of the Area of Integration 

Since S{x) is a non-negative function, the asymptotic expansion of 
the deformed heat trace 'c)q>{t) will be determined by a neighborhood of points 
where S{x) = 0, i.e. by M^: 

Lemma 3.1. There holds as t 0^ 

e$(t) ~ {inty^ j dx exp . 

Me 

Proof. Let p, i?, C be such constants that the estimation ()3.2|) holds. Let Br 
be the ball with radius R centered at the origin: Br = < x : cr (a;, 0) ^ — >. 
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We split 0$(t) into the sum of three integrals: 

e^{t) = hit) + hit) + hit) , 



where 



hit) =(47rt)-i j dx exp (-^^(a^)) , 

/..e.p(-i..), 
hit) =i^nt)-' j rfx exp . 



Me 



hit) =(47rt) 

Br\M, 



We will show that hit) and hit) are exponentially small as t — 0"*^, i.e. they are 
asymptotically equivalent to zero and the asymptotics of is determined 

by hit) only. 

Since \ Mi; is compact and Six) > on this set, it is separated 
from zero: Six) ^ 5 > for all x G \ M^. Hence, as t — 0+, we have 

hit) ^ (47rt)-^ TxR^ exp ( -— ) ~ . 



2t^ 

For /3(t) we use (j3.2j) to estimate the integrand at infinity: 
/3(t) ^(47rt)-i j dxexp(^-ic[a(x,0)f^ 



5fl 
OO 



= (2t) ^ j rdr exp 



Cr^p 
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where T{x, a) is the incomplete r-function: 

oo 
a 

It is well known (see, e.g. [H]) that as a ^ oo 



r(x,a) = a"-^e-" ( 1 + ( - ) ) . (3.6) 



a 



Hence, we have as t — 0^ 

Therefore, e$(t) ~ /i(t). □ 
So, after the change of coordinates = x"{y, z) in we get: 
e<,(t) ~ (47rt)-i y dvol(y) y t/^ exp (^-^^S{x{y, z))^ (1 + z x{y)) . 

3.1.3 Taylor Series Expansions 

By scaling the variable z — *• y/tz we obtain 

e/Vt 

Q^it) ~ y" dvol(|/) y dz exp (^-^5 (x (y, ^tz^j^ {l + v^zx(2/)) . 

Now we expand the integrand in the last formula in the power series in \/t. 
Below we compute the series expansion for S* (x (y, v^2;)). 

We define 

^"(a;) = (T"($(a;), x) = $"(a;) - x". 
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Then 



Since 2; = is the fixed point set of the mapping $, we have 



^"(/(2/)) = 



Therefore, 



k=l 



where 



gk 



2=0 



Qk^a 



dx"^ . . . dx"^ 



2=0 



Now we obtain 



S [x [y, Vtz^^ ~ IT^*^*^^'* 

fc=2 



where 



1=1 ^ ^ 



■I)- 



We consider a non-degenerate case (see the assumption ($.3)), when 



2 = 



7^0 



for any point of S. In this case 



S2 = > 



everywhere on E. 
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Remark. It is worth mentioning, that if S2 = somewhere on E, that is, 
the first non-vanishing term of the Taylor series ()3.7|) of S{x{y, z)) is of degree 
higher than 2, then the asymptotics of Q<s,{t) as t ^ 0"*" will be substan- 
tially different. Note also, that since S is compact, this restriction means that 
S2 ^ (5 > for some constant 6. 

The next several coefficients are: 

53 = 3^fi)^,(2), (3.9) 

54 = 4^fl)^„(3)+3^^2)*„(2), (3.10) 

55 = 5 *fi)^„(4) + 10 ^f2)*a(3) , (3.11) 
Se = 6 *fi)^„(5) + 15 ^f2)^a(4) + 10 *"3)*a(3) • (3-12) 

Now we rewrite the integral for 0#(t) in the form 

e/Vi 

e$(t) ~ ^ y" dvol(|/) j dzexp(^-'^z'^B{t,y,z), (3.13) 

where 

/ I °° +k/2-l \ 

B{t, y, z) = exp f -- J] -^s, (l + t^'^zn) . (3.14) 
Our next step is to expand B in the power series in 

^1/2 

oo oo 

^-^ftfct'^ + t^/'^M^ (3.15) 

fc=0 fc=0 

and compute integrals over z. It is easy to see that in this Taylor series: 
1. All coefficients hk and hk are polynomials in 2, x, and s^, ^ 3. 
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2. Integer powers of t have coefficients bk with even powers of z only. 

3. Half-integer powers of t have coefficients with odd powers of z only. 

All bk will disappear after integration over z, since they are odd functions, so 
we do not need to keep track of them. For computing integrals involving bk we 
will use the following result. 

Lemma 3.2. Let a > 0, G Z, ^ 0. Then there holds as t — O"*" 

e/Vi oo 

dz e-'^^'z^'' I dze^^'^z' 



Proof. We have as t — 0''" 

oo e/v^ 

jdze-'^^'z''- I dze^^^'z'' = a~^''+'y'r 1^^^,'-^^ r^O. 

-°o -e/Vt 

Here we used the asymptotic expansion ()3.6p . □ 

Therefore, by substituting the expansion ()3.15|) into ()3.13j) and using 
Lemma (3. 2t we obtain 

oo 



k=0 



where 



Ak= dvol(y) ttk 



E 



and 



4vr 



oo 

^ J dz exp (^-^ z^j 6fc . (3.16) 
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Due to Lemma 11.11 integration over z will bring inverse powers of S2- This 
completes the proof of Theorem Kill □ 

Now we will compute the first several coefficients Ofc. By expanding 
the function B{t, y, z), given by (j3.14p . in the power series in t^/^, we obtain 

&o = l, 

^^='-12^^"-48'V' ' 

1 1 \ 6 1 1 1 2^ 8 

240 1440 / V576 2880 4608 V 



+ 



-S3X 



-S3S4 ) z^^ 



4^12 



10368 13824 ^ V 497664 

Therefore, by using formula ()3.16p . Lemma ITTTl and performing straightforward 
computations, we get the following lemma. 

Lemma 3.3. Let conditions of Theorem V^ . 1\ he satisfied. Then 
ao = (47rs2)~^/^ , 



1 

4' 



02 = (47rs2)-^/' 
35 



1 ^ 
12' 



35 



3 2 



2 5 
^2- + 


3'='2 




^35 


b) ^^2-^ + ( 


^12 


385 4 6 




288'3'2 





7 ^ 
12' 



35 
96' 



4 / ■^2 



where coefficients Sk are given by (j3.8|) - ()3.12p . 



3.2 Zero-Dimensional Submanifolds 

In this section we consider the case when the fixed point set S is just 
a point, that is, dimS = 0. 
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Theorem 3.2. Let ^ be a smooth mapping $ : ^ M^. Let S be its fixed 
point set: S = {x G : = x}. Let S be a real-valued function 5 : — > M 
defined by S{x) = (t($(x),x). Let $ satisfy the assumptions: 

($.1) H is a one-point set, S = {.^o}- 

($.2) There exist constants p > 0, R > 0, C > such that for any x G , if 

a (x, 0) > — , then 
2 

S{x) ^C[a(x,0)f . 
($.3) .^0 non- degenerate critical point of the function S. 
Then there exists the asymptotic expansion as t 

oo 

e^(t)~ J^t'^Afc, (3.17) 

k=0 

where scalars depend polynomially on derivatives of the function S and the 
inverse of its Hessian matrix, evaluated at C,o- 

Before the proof of the theorem, let us note that, given the fact that 
coefficients Ak are polynomials in derivatives of the function S, it is possible to 
write down their general form by using dimensional considerations jTTj. The 
deformed heat trace is dimensionless, so the expansion on the right hand 

side of (j3.17p must also be dimensionless. From the heat equation, t has the 
dimension L^, so the dimension of A^ is The function S has the dimension 

of (L — dimension of length), since it is equal to the world function. Its 
m-th derivative with respect to the space variables has the dimension L^"™. 
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In particular, its second derivative is dimensionless, and all derivatives of order 
higher than two have dimensions of negative powers of L. 

In order to construct scalars A^, it is necessary to contract all indices 
of derivatives with indices of the inverse of the Hessian matrix of the function 
5, which we will denote by Q. The order of this contraction can be different, 
therefore, the coefficients Ak are linear combinations of invariants of the form 

T(m, V) = (g^.-.^g)"-"" ® . . . ® V-5)^^^^^_^^^^^^ , 

M 

where m = (mi, . . . ,m;) is a multi-index, |m| = mi + . . . + = 2M, and 
^ S2M is a permutation of numbers 1, . . . , 2M. Note, that mj ^ 3 for any i, 
since Aj^ depends polynomially only on derivatives of the function S of order 
higher than two (besides, the function S itself and its first derivative vanish on 

Let us introduce the length, the order, and the dimension of the term 
T (we omit arguments): 

len T = I , 
ordT = 2M, 
dimT = 2M-2/. 

In other words, the length is the number of derivatives of the function S in the 
term T, the order is the sum of orders of all derivatives, and the dimension is 
the (negative) space-dimension of T. Clearly, 

dim T = ord T — 2 len T . 
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Now we are ready to write down the general formula for the coefficients A^: 



2k 



Ak = Y^ c{m,(f)T{m,(f) 



1 = 1 m={»ni,...,»n,) (peS2(k + l) 
|m( = 2(fc + i) 



where c(m, if) are some numerical coefficients. Note, that this sum is finite for 
any k. 

Proof of Theorem VJ.^A The idea of this proof is very similar to the one of The- 
orem 13.11 but in this case we need to use the multi-dimensional version of the 
Laplace method. In the same way, due to Lemma 13.11 we can reduce the area 
of integration to a small neighborhood of S and introduce there a system 
of Cartesian coordinates with the origin at .^o- Then we expand the function 
S in the Taylor series in scale the variables z" — \/tz°' and obtain 

e.WM4.)-'/<<^.exp(-is„,.".^)s,M). (3.18) 

where 

^p5.,....,^^^..z^O, (3.19) 

s\ . 

Again, B has the expansion 

oo oo 

B^Y.^ut'' + t^'^Y}kt'' (3.20) 

k=0 k=0 

with coefficients and being polynomials in and derivatives of S of order 
higher than two. And there is no need to keep track of b^ since they do not 
contribute to the asymptotic expansion of the deformed heat trace 0<i,(t). 



Bit, z) 


= exp ^ 


oo 
k=3 




d 


d 


dz"^ 


dz"" 
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Since is a non-degenerate critical point of the function S, its Hessian 
matrix Saj3 is invertible. We will denote the inverse matrix by Q. It satisfies 
the equation 

By substituting ()3.20|) into the integral ()3.18|) and applying Lemma 11.21 we 
conclude that 



k=0 



where 



Afc = ^ y d'z exp S^pz'^z^^ h . (3.21) 

Due to Lemma fl.2l integration over z will bring factors of Q. This completes 
the proof of Theorem 13.21 □ 

Now we will compute the first several coefficients A^. In order to sim- 
plify formulas, we will use the following compact notation for the (symmetrized) 
derivatives of the function S: 

S{k) = Sf^y^,„y^) , 

~ k\ '"i---'"* ■ ■ ■ 
By expanding the function B{t, z), given by ()3.19p . in the power series in t^/^, 
we obtain 

&o = l, 

b2 = -2%) + 4 %)'5(3) + g S{A)S(^4) - — 5(4)5(3)5(3) + — 5(3)5(3)5(3)5(3) . 
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In terms of the mapping $ and its derivatives, derivatives of the functions S 
are 

Sial3) =gaf3-2 9^,(a^^p) + , (3.22) 

= -3 ^M(«'^''/37) + 3 g^.u^^'ia^'^dj) , (3.23) 

Sia^^S) = -4 ^^(a$^/375) + 4 Q f.u'^^ (a'^'' pyS) + 3 9 ^,u'^^\ap'^%S) , (3.24) 

'5{a/375s) = -5 gf.ia'^^p^Se) + 5 9 [a'^" p^Ss) + 10 ^^^$^(a;3$%fc) , (3.25) 
S{aP'y5e(;) = 9 ^l{a^^ p-ySeQ + ^ 9 /lu^^ {a^'^ p^ySeQ + 15 5'^i,$''(q/3$''^5£^) 

+ 10(7^,<I>^(,^^<I>V)- (3-26) 



These formulas can be obtained as a particular case of (j4.16p - (j4.2(jp . where all 
derivatives of the world function vanish except for the second order. 

Let T : Vm+2 — ^ Vm be the operator of contraction of a symmetric 
tensor with the matrix Q: 



al3ui..M„ 



By using formula ()3.21|) . Lemma ll.2t and performing straightforward compu- 
tations, we get 



Ao = VdetQ 
Ai = v^det Q 

A2 = ^/detQ 



(%) + (^(3) V ^(5)) + I V %)) 



35 



385 

288 



(S'(3) V S(3) V 5(3) V 5(3)) 
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Since expressions involving symmetrization are inconvenient for practical com- 
putations, we will expand them in a symmetrization-free form. We will use the 
following notation 



(3.27) 

that is, the number in brackets is the number of indices contracted with Q. By 
using corollaries 12.11 12.21 12.31 12.71 and 12.61 we obtain the following lemma. 

Lemma 3.4. Let conditions of Theorem VJ.^A he satisfied. Then 

= Vdet Q , 

Ai = ^/detQ 



A2 = Vdetg 



1 1 ^^^^ 1 



I ^ c c /n VI !^5^ 1/21^6/0 ^3 ^7 /nail's ''" Qr , Q, , Q, ^D'^i'^a 

I —'->iyiU2U3iy4'Ju5U(iL'7l's'^ ~'J[2]l/i'J[2]l/2*-'[4] V 

— -;^S[2\uiSy^i,.^y^S[2\u.^yQQ^^^^Q^''^^''Q^*^^ — -^S[2\uiS[2\u2S[2\u3UiQ'^^^''^Q^ 

_ 1 Q, , Q Q r^i'2v% r^yiyi r^yiys, 

2 'J[2]vi*-'l/2!^3!^4'-^I^5!^6!^7l^8V V V V 

C 0^^/01^11^4/01^21^5/01^31^6 

2^^I^lI/2l^3'-^!^4l^5l^6'-^[4J V V V 

_ 1 c Q , r^yiyif^y2ybr^f3fTr^ffifs 

_ ^ c c c /ni^ii^4/ni^2i^7/ni^3i^8/n!^5i^9/n!^6i^io 
2 '-'!^i!^2i'3 '-^!^4!^5i'6 '-^'^Ti^si'gi^io 

+ ^S[2\u^S[2\u2S [2] 1/3 5'[2] 1/4 Q Q 
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24 V V V 

+ Si2] ui Sy^ ^3 S\2\ 1/5 5'[2] Q^^^'^Q ^■^ '^^ Q 

~r ^^[2\Ui^U2l'3U4^Ul,h>(iU7^Ush'9l^l0^ ^ ^ ^ ^ 

-^2 ^^'^^''■^ ■''^^'^^ ^^"'^^^'^^•^ V V V V V 

I ^ c c c c /ni^i»^4/ni^2»^5/ni^3'^7/ni^6'^io/ni^8i^ii/ni'9i^i2 
~r ^ '-'1^4(^51^6 '-^i^7'^8'^9 '^io'^ii'^i2 ^ >ct; 

i^c c c c /n!^i!^4/n!^2!^7/n!^3i^io/ni^5!^8/ni^6'^ii/ni'9!^i2 
o 

where tensors S[2k]ui...u,n ^-re given by (j3.27|) . (j3.22|) - (j3.26p . and Q is the inverse 
matrix to Sap- 



CHAPTER 4 



Spectral Geometry of Submanifolds: Curved Manifold 

In this chapter we will consider a curved manifold M, and a smooth 
mapping <I> : M — M on it, such that its fixed point set S is zero-dimensional, 
that is, it is just a single point. We provide a new proof for asymptotic expan- 
sion Theorem 14.11 This proof will give us a way to compute the asymptotic 
expansion coefficients. The coefficients Ai, in Lemma and the coefficient 
A2 in Lemma [4.31 are computed explicitly for the first time. 

4.1 Main Theorem 

Theorem 4.1. Let {M,g) be an n- dimensional compact smooth Riemannian 
manifold without boundary. Let ^ be a smooth mapping $ : M — » M. Let S 
be its fixed point set: S = {a; G M : = x}. Let $ satisfy the assumptions: 

($.1) is a one-point set, S = {xq}. 



k=o 

where the coefficients A^ are scalar invariants depending polynomially on co- 
variant derivatives of the curvature of the metric g, symmetrized covariant 



($.2) det(/-rf$)(xo) ^ 0. 



Then there exists the asymptotic expansion as t — > 0"*" 



00 




(4.1) 
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derivatives of the differential of the mapping $, and the matrix 

Q = {{I-d^fg{I-d<^)}-\ (4.2) 

evaluated at xq. 

Proof. The idea of this proof is the same as for a flat manifold. We will choose 
a convenient system of coordinates, reduce the integral over the whole manifold 
M to the integral over a neighborhood of the fixed point set S, expand integrand 
in Taylor series at Xq, and compute the asymptotics of the integral by reducing 
it to the standard Gaussian integrals, which are known exactly. 

In this proof we will make several references to some computations in 
the upcoming sections. However, it does not affect rigor, since those computa- 
tions do not use this theorem. 

Let 5 be a real-valued function, defined in some neighborhood of S 
by S{x) = a (x, Obviously, Xq is a point of minimum of the function 

so S{xo) = and S.^ixo) = 0. We will show below (Section 1021 (ElSl)) that 

SM^o) = [S^a - {d^xo)ra}g,A^o){S''f, " (d^xo))" p] , 

or, in the matrix notation, 

VVS{xo) = {I- d$ (xo)}^^?(xo){/ - rf$(xo)} . (4.3) 

Therefore, the matrix Q defined by (|4.2p is nothing but {VVS'(xo)} ^■ 

From now on for simplicity we will omit a semicolon in covariant 
derivatives of scalar functions, for example, Sa = S-a, Sap = S-ap, etc. 
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Let Mg be an open geodesic ball with center at the point Xq and small 
enough radius e, so that the functions a and S are well-defined in M^. Since 
M \ Mf, is compact, there exists a positive constant 5, such that for any point 
X G M \ we have 

dist(x, $(x)) > 5. 

From the last inequality and the standard off-diagonal estimates for the heat 
kernel in ^2] it follows that as t — >^ O"*" 

dvol(x)f/(t;x,$(a;)) ~ 0. 



Therefore, 



e$(t) ~ / dvol(x)f/(t;x,$(x)) 



and according to ()1.8|) we have 
e$(t) ~ j dvol(a;)(47rt)-"/Mi/2^x, $(x)) exp ("^^ ^(i; $(a;)) . (4.4) 

Me 

Now we will choose a convenient coordinate system in and expand 
the integrand of ()4.4j) in the covariant Taylor series at Xq. 

4.1.1 Normal Coordinates 

Let X be a point in M^. Due to ()1.6|) . the components of the tangent 
d 

vector Z = to the geodesic connecting the points Xq and x at the point 

xq in the direction of the point x and with the norm equal to the length of this 



geodesic are given by 



OXq 
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According to (|1.7p . we have 

X = exp^„(Z) . 

One can use this vector to define new coordinates in by assigning to a 
point X with coordinates new coordinates z°' equal to the components of 
the vector Z in the coordinates x'^. Let (f : R" be such a coordinate 

mapping. Since is a geodesic ball with radius e, the image of is 

(p{M,) = {zeW : \\z\\ < e} . 
The differentials are related by 



dz- = ^dx^ = -,-(xo)^^(xo,x)cix^ 



Thus, for the measures we obtain 

- g^/3g^7 (^o,a;)J dx 
= \g{xo)\~'^^^ A{xo,x) \g{x)\^^'^dx 
= \9{xo)\~^^^ A{xo, x) dvol(x) . 



So, the Riemannian volume element has in new coordinates the form 

dvol(x) = A~^{xo,x) \g{xo)\^^'^ dz . 

After the change of variables x —>■ z integral ()4.4|) takes the form 



|^Xoj|^''^exp 



S{x{z)) 



\\z\\<e 



2t 

X A-^{xo,x{z))A^^^{x{z),^x{z)))n{t;x{z),^x{z))) . (4.5) 
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4.1.2 Taylor Series Expansion 

Let / be a smooth function on the manifold M. Then in the neigh- 
borhood Me of the fixed point xq this function can be expanded in covariant 
Taylor series jS] 

oo ^ 

/(^(^)) = E^(V'/)(.......)(^o)^-...^-, (4.6) 

k=0 

where are the coordinates introduced in the previous section. Since this ex- 
pansion involves only covariant derivatives of the function / and the coordinates 
are covariant derivatives of the world function, all parts of this expression 
are written in an invariant form. Note also, that all covariant derivatives of the 
function / are symmetrized, which makes (|4.6p very convenient for practical 
computations. We will use the notation 

/(fc) = Sym (VV) (xo) , 

We will expand each factor of the integrand in ()4.5jl in the Taylor 
series at Xq separately and then multiply the expansions. Note, that we will 
need to compute symmetrized covariant derivatives of functions of the form 
/(x, $(x)) at the point Xq, which has the property $(xo) = Xq. 

The function Vt has asymptotic expansion ()1.10|) in powers of t. Let 

oo 

a^{x{z)) = u;rr,{x{z), ^{x{z))) = ^am{k) , 

A;=0 

that is, the first index stands for the coefficient of the power series in t, and the 
second one for the coefficient of the covariant Taylor series in 2;. It is convenient 
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to represent the Van Fleck-Morette determinant in the form (2] 

A = exp(2C) . (4.7) 

The function has a known expansion in a Taylor series with the first two 
coefficients being equal to zero 0, 

C,{xo,x{z)) = J^C(fc) • 

k=2 

Let 

V{x{z)) = C(x(z),$(x(z))) 

and 

oo 
k=0 

We will show in Section 14.2.51 that r/(o) = ^(i) = 0. The Taylor series of the 
function S also starts from the second term, 

oo 
k=2 

4.1.3 Scaling of Coordinates 

Now we rescale the coordinates 

^ Viz"" . 

Then all terms in covariant Taylor series are transformed as 



f{k) ^^^^/{fc) , 
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and, hence, they become asymptotic expansions in t^/^. By substituting ex- 
pansions from the previous section and p.lO|) into ()4.5p . we obtain 

e^{t) ~ (47r)""/2 f dz |^7(xo)r/'exp (-^ S^^)) B{t) , (4.8) 



\z\\<£/^/t 



where 



ml 

k=l / m=0 k=0 



and 

^fc = ^ S(k+2) + 2 C{fc) - Vik) ■ (4.9) 

Note, that terms S(^k) and ?7(fc) depend on the mapping $, while does not. 
Let fefc be the coefficients of a power series in t^/^ for the function B, 

oo 

B{t) = Y,t^'Hk. (4.10) 

fc=0 

It is easy to see, that the coefficients hk are polynomials in z. Moreover, they 
are polynomials in hk and 'dm{k) ■ From |3] and our computations in Section 14.21 
it follows, that and am{k) can be expressed in terms of covariant derivatives 
of the curvature of the metric g and symmetrized covariant derivatives of the 
differential of the mapping $. Note, also, that half-integer (integer) powers of 
t correspond to polynomials of odd (even) degree in z. By substituting ()4.10|) 
into ()4.8j) . we obtain 

e^it) ~ Y.^'/' (4^)""^' / l^?(^o)r/'exp (^-i 5(2)) 

||z||<e/v^ 



(2) . 
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In the same way, as in Lemma 13.21 we expand the integration over the set 
||2;|| < el\/t to the integration over the whole M", i.e. one can show that 



y dz exp 1^-^ S'(2)^ z^^ ... z''^- ~ y exp 1^-^ S'(2)^ z'K . . £ 



k 



\\z\\<e/Vt 

for both odd and even k. After this operation all terms with half-integer powers 
of t disappear, since they contain integrals of odd polynomials over a symmetric 
set. Therefore, we have 

e$(t) ~ f^t' (An)--/' I dz |^?(xo)r/'exp (-^^(2)^ 62,. 



By comparing this expression with (j4.H) . we finally obtain 

Ak = (47r)-"/2 J dz |(7(xo)r/'exp (-^^{2)) hu ■ (4.11) 

In this expression: 

1. S'(2) = -^Sapz'^z^ , where Sap is the positive-definite quadratic form, de- 
fined by 

2. &2/c are polynomials in z, defined by (|4.10|) . 

3. All integrals in ()4.1H) are standard Gaussian integrals which can be com- 
puted by using Lemma 11.21 

4. are polynomials in covariant derivatives of the curvature of the met- 
ric g and symmetrized covariant derivatives of the differential d^ of the 
mapping $. Due to Lemma II. 2| Gaussian integrals in 1)4.111) depend on 
the inverse of the matrix S^p, that is, the matrix Q. 

This completes the proof of the theorem. □ 
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4.2 Computation of the Coefficients Aq, Ai, and A2 

In the remaining part of the chapter we will develop an algorithm 
for computation of the coefficients Aj. in terms of covariant derivatives of the 
differential (i$ and the curvature of the metric g. We will compute coefficients 
Aq, Ai, and A2 explicitly. In the next section we describe tools for construction 
of Taylor series of two-point functions. 

4.2.1 Coincidence Limits of Two-Point Functions 

Let f : M X M ^ R he a two-point function defined on a manifold 
M. We will denote by [/] its coincidence limit, that is, 

IfK^) = lim fix,x') . 

Obviously, covariant derivatives of the function / at different points commute 
with each other. We will use primes for indices of covariant derivatives at the 
point x'. For coincidence limits of mixed derivatives of the function / there 
holds [m 

where dots stand for a fixed set of some indices and may include other covariant 
derivatives. By using this property repeatedly, we obtain 

[f..,a'P'] = [f...U - [f.-M - + {f...;Po] , (4.12) 

[/...;a'/3'7'] = [/...];a/37 ~ [f...n];af3 — [f...;f3];a-f " [f...;a\;(3-y (4-13) 

..;7a];/3 "I" [/...;/3q];7 [/...;7/3a] • 
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One can easily obtain the following general formula of this kind for [/...;^'^...^^]. 

Let Nm = {1, . . . ,m} and P(Nm) be its power set. For any I G P(Nm), such 

that / = {zi, . . . ,ik} and ii < . . . < ik, let be the increasing sequence 

> 

of indices ■ ■ ■ Hi^. and be the decreasing sequence of indices /ij^. . . . Hi^. 

i 

Then 

where J = N^, \ /• 

4.2.2 Derivatives of the Deformed Diagonal Functions 

Let / be a two-point function / : M x M ^ M. In this section we 
will develop an algorithm for computing symmetrized covariant derivatives of 
the function g{x) = /(x, $(x)). We will agree, that indices with primes denote 
covariant differentiation of the function / with respect to the second argument 
and indices without primes — with respect to the first argument. Obviously, 
derivatives with respect to different arguments commute with each other. 

Let be the set of two-point tensors of the type (p, q) at a point 
X and of the type (r, s) at a point x', that is, the set of linear mappings 

T : T^M X ... X T^M x T*M x . . . x T*M 

9 p 

X T,/M X ... X T^,M X T*,M x . . . x T*,M ^ R . 
' V ' ^ ' 

s r 

If some of the numbers p, q, r, s are equal to zero, we will omit them. The 
external tensor product of two-point tensors A G and B G ^^W^ is the 
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two-point tensor AMB e gXl^l+l^ defined by 

[/i ^ ^J/3l.../3,,.i...i/„5i...5^7rJ...< - ^0i...l3,8[...5'^ ^ui...Uu7T[...n'^ 

We define a mapping • : ^W^ x ^VV ^ p+r^V by 



It is easy to see, that the n-th symmetrized derivative of the function 
is a hnear combination of terms of the form 



T{k,m;j,, ...,3k) = Sym (V'^'V"/) • ((V^M$) 
where numbers k,m,ji, . . . ,jk are non- negative integers, such that 

k + m + ji + ...+jk=n. (4.14) 

A simple computation shows that there exist 2" such terms. 

Let numbers c{k,m;ji, . . . , j^) be the coefficients of this hnear com- 
bination, that is. 



n n—k 



m; ji, . . . , jfc) T{k,m;ji, . . .,jk) ■ 

k=0 m=0 j=(ii,..-jfe) 

\j\—n — 7n — k 

We will derive a recursive formula for these coefficients. Note, that the term 
T{k,m] ji, . . . , jk) of the n-th derivative can be obtained from the following 
terms of the (n — l)-th derivative: 

1. T{k,m — 1; ji, . . . ,jk), by differentiating / with respect to the first argu- 
ment, if m > 1. 



,Jk) 
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2. T{k — l,m; ji, . . . ,jk-i), by differentiating / with respect to the second 
argument, ii k ^ 1 and jk = 0. 

3. T{k,m;ji, . . . , - 1, . . .,jk), by differentiating ci$, if jj ^ 1. 

Let c{k,m;ji, ■ . . ,jk) = 0, if one or more arguments are negative. The above 
rules give us 

k 

c(/c,m; ji, . . . , jfc) = c{k,m- 1; ji, . . . Jk) + ^c(/c,m; ji, . . . , ji - 1, . . . 
if jfc > 0, and 

k 

c{k,m;ji, . . . , jfc) = c{k,m- 1; ji, . . . , jfc) + ^c(/c,m; ji, . . . , ji - 1, . . . 
+ c(A; - l,m; ji, . . . , j^-i) , 

if jk = 0. The obvious initial condition for this recursion is 

c(0,0;) = l. 

One can also get a non-recursive formula for the coefficients c{k,m;ji, . . . ,jk)- 

Lemma 4.1. Let n G N and k,m,ji, . . . ,jk be non-negative integers, satisfy- 
ing (jUIlD. Then 

qft,m, ji, . . . ,jk) _ _ . + j.)(2 + J, + ...{k + jk + ...+3,)- 

Proof. It is not hard to check that the non-symmetrized covariant derivative 
9ui...u„ is the sum of all terms of the form 

/ i^,p(i)...i',p(„) 
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for all possible values of fc, m, ji, . . . , and all permutations of numbers 
l,...,n satisfying certain conditions, that we will determine below. Sym- 
metrization of indices z/i . . . z/„ makes all terms with fixed k,m,ji, ■ ■ ■ ,jk the 
same. Hence, the coefficient c{k,m; ji, . . . , jk) is the number of all possible 
permutations ip. 

Let Jj be defined for z = 1, . . . , /c by 

Ji = m + i+ji + ...+ ji-i . 

Note that indices . . . i^^(m) correspond to the derivative of /, index z^i^(ji) 
corresponds to the z-th factor d^, and indices . . . correspond 

to the derivative of this factor. This observation leads us to the following 
conditions on the permutation ip: 

1. (/?(!) < (p{2) < ... < ip{m) ; 

2. <p{Ji) < <p{Ji + !)<...< ^{J, + ji) ; 

3. <p{Ji) < (^(Ja) < . . . < <p{Jk) . 

In words, indices on every factor and the first indices of factors must be 
in the "increasing alphabetical order." So, there are (^) ways to choose in- 
dices corresponding to the derivative of /. After this step </?(Ji) is determined 
uniquely and there are ("~™^~^) ways to choose indices corresponding to the 
derivative of the first factor d^. By continuing this reasoning, we get 




X 



( 



n — m — 2 — ji 



32 



) 



( 



n - m - {k - I) - ii - . . . - jk-2 



Jk-i 



) 
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After straightforward transformations, we get the desired formula. □ 

4.2.3 Symmetrized Derivatives of the Function S 

Now we will apply the technique developed in Section ll.2.2l to the com- 
putation of symmetrized covariant derivatives of the function S{x) = a{x, 
up to the sixth order. Since we are interested in values of the derivatives on 
the fixed point set S, all appearing derivatives of the world function will be ac- 
tually equal to their coincidence limits. We will take into account the following 
properties of the world function 0: 

1. The world function is symmetric in its arguments. 

2. The first two derivatives with respect to the same argument commute. 

3. [a] = 0. 

4. [a^] = 0. 

5. [O-Q/?] = -[(Tal3'] = QaP- 

6. [0"q/37] = = 0. 

7. Kmi...^.,)] = 0, ^ 3. 

8. [(TQ(^.i...f.fc)] = = 0, A; ^ 2. 

The recursive algorithm, described above, for computing the coef- 
ficients c{k,m; ji, . . . , jk) was realized as the Python script listed in Appen- 
dix lA.ll This script also generates preliminary expressions for derivatives of 
the function S in the I^TgX format. 
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For the derivatives of the differential of the mapping $, evaluated 
at the point Xq, the notation 

will be used. The first six symmetrized derivatives of S are 
5 = 0, 

Sa = 0, 

S(al3) = 9al3 " 2 (a/S) + '^n{a'^'^l3) , 

S{al3-y5e) = "5 ^'(a/3'y5e) + 5 /^{a'^'^ p-ySe) + 10 (^{als'^'^'ySe) 
+ 30 [a^V'(a/3]^''7*"5£) 

+ K'u'^'ia]{5 ^^fS^%<il^Se) + 25 ^'^^^^^^^.J } 
+ 10 + 10 [a^^VC^] *''7*"5*^e) , 

S{aP'ySeO ^ ^(a/37<5eC) + 6 ^'^{a^I^''/375eC) + 15 ^ ij.{ap'^'^ ^SeQ + 10 ^^(Q/37^''<5eC) 
+ [a^V'MljeO *^^''5,C) + 45 ^%5^"eC)} 

+ 54 ^^^^%fe^^c) + 48 ^^;37^'^fe^^C) } 

+ 18^^(„^^^^^«fe^\) + 15^^(„^^%^«,^\)} 

+ 60 [a^V'(a/37]^''5^''eC) 
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+ 42 ^^^^%*«,*\)} 
+ 15 [a^ 

The second derivative can also be written in the form 

= (5^ - {d<!>ro)g,u{S''p - {d<^>y^) . (4.15) 
We need the following coincidence limits for derivatives of the world 

function: 

[o'lj.'u'iafSyS)] , [0"(AjV'f )(a/37)] • 

We will reduce them to the limits which are known j3] 



o 7 



For this purpose we will use properties of coincidence limits (I4.12|) . ()4.18j) . 
expression for the commutator of covariant derivatives (jl.4j) . symmetry prop- 
erties of the curvature tensor and Bianci identities (|1.2|) . (jl.3|) . After 
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straightforward but somewhat cumbersome computation we have: 









W ^JL' u' a] = 






[o'^i'u'^'tt'] = 






[o'n'u'Ca/s] = 


1 


5 1 

- R{fi\a\u)l3;i ~ 12 ' 


[o'^'u'CTT'a] = 


5 

g R(p\a\u)(t,7r) - 


- g Rii\a\n){^,■,u) , 


[o'l^i'u'iafS'rS)] = 


2 

--^R^l{a\ly\(^nS) 


g 

- — Ra{a\^l\|3R''-/\u\S) , 




(-^R^\a',; 


^) , _ 1 o eokl t; 

7) 3 (a p 



("ICI/3){m;!-) 



Detailed computation is given in Appendix IA.31 

By combining obtained expressions, we finally have 

SiaP) =9ap-2 -^^^13) + -^f^ia^^p) , (4.16) 
^(a/37) = -3 *(a/37) + 3 ^,^{a^''|3f) , (4.17) 

(4.18) 

+ 10 ^^.iaf^^^-ySs) - 5 ^'(a/375e) + 5 ^ p^a"^^ f^^Se) , (4.19) 



-^9i?„ 



3 



Rfj.{a\u\l3;\^TT\ + 8 Rg^u{aR^\^TT\P \'^'^ -y"^" s'^K'^'' 
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+ - {6%|^|,.)(a;/3 - 10 R^,{a\u\m " ^?(a|M|/3;|i^| " ^5(a|i.|/3;|M| } 

+ 9 ^^^^ V*^C) + 8 *V*'5e*^c)} 

(4.20) 

4.2.4 The Heat Kernel Coefficients 

By using the general algorithm for computing coincidence limits of the 
heat kernel coefficients and their symmetrized covariant derivatives, developed 
in 13 El, we obtain 

[^o] 

[l^l;a/3] 
[^2] 



6 



R. 



12 



R-a 



1 



36 



R^ 



— R/ 
15 



1 

90 



90 
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where 

(4.21) 

So, for the coefficients afc(x) = ujk{x, at the point Xq we have 

ao(a;o) = 1 , 

ai{xo) = -\r, 
o 

ai-a{xo) = -— - — R-^'^^a , 

4.2.5 The Van Fleck-Morette Determinant 

For the function ( in representation ()4.7|) we have [3] 

[C] = o, 

[Ca]=0, 
D 

[C{a(3-y)] = ^ -R(«/3;7) > 

3 1 

In the same way, as for the world function, we obtain non- symmetrized deriva- 
tives of the function ( from these expressions for symmetrized derivatives. After 
performing straightforward transformations, we have 

[(0/37] = Y -^(a/3;7) ' 



[C 
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where 

3 1 1 

FaMS = -^(<^/3;7'5) + YS ^) ~ 9 -^M(a-R^/3|5|7) • (4.22) 

By using the method developed in Section l4.2.2| for derivatives of the 
function ri{x) = ({x, $(x)) at the point Xq we obtain 

r]{xo) = , 
r]a{xo) = , 

11 1 1 

2 1 r 1 

+ R{aP;\iM\'^'^^5) + '2 Rij.u;{a'^'^l3'^''^5) + -R/i!/;(a/3^^7^'^<5) 

— 2 -R(«/3;|m*^75) ~ Rfiia;l3'^'^-y5) — R^lu■,{a^^ P^"" -yS) 

3 

3 1 
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4.2.6 CoefRcients bk 

By direct expansion we get the coefficients Bq, b2, and 64, defined 
by ()4.10|) (all objects are computed at the point Xq and we will omit arguments) 

&o = l, 

h = -am - h + ^hl, 

64 = 2 ^2(0) — «i{2) + + — - h\) ai(o) 

11 1 

- /i4 + /il/i3 + 2 ^2 - 2 ^1^2 + ^ ^1 • 

Due to ()4.9p and taking into account that C(i) = ^(1) = 0, we have 

hi = -5(3) , 

^2 = - <S'{4) + 2 C(2) - ^(2) , 
^3=2 -^(S) + 2 C(3) - ^(3) , 
^4=2 -^(6) + 2 C(4) - ^(4) • 

Therefore, 
fco = l, 

1 ~ 1 /- \2 

62 = -ai-2 C(2) + r/{2) - 2 %) + g (^"^Cs) ) , 

^4 = - 02 — «i{2) + 2 aiC(2) — ai''7{2) + aiS[4) — 2 C{4) + Vii) + 2 '^'(3)01(1) 

/~\2 1_ 2 ~_ 1~ 1~~ 1~_ 

+ 2 (^C(2) j + - (^(2)) - C(2)^{2) - - 5(6) + 2 %C(2) - ^ 5(4)7/(2) 

1 /- N 2 _ _ 1 ~ _ 1 ~ 1/^x2 

- g Ol (^5(3) j + 5(3) C(3) - 2 '5(3)^(3) + ^ 5(5)5(3) + - (^5(4) j 

(^(3)) C(2) + ^ (5(3)) ^(2) - ^ 5(4) (5(3)) +^ (5(3)) • 
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4.2.7 Gaussian Average 

Since integrals in ()4.11|) are Gaussian, it is convenient to introduce a 
Gaussian average for a function /, defined on M", by 

(/) = (4vr)-"/2f,-i I dz |^7(xo)r/'exp (^-1 S^^z'^z^^ f{z) , 

where 

Q = |det (/ - d^{xo))\~'^ . 

The Gaussian average is normahzed so that (1) = 1. Recall that -Sa/sz'^z^ = 
S(2), SO ()4.11|) takes the form 

Ak = g{b2k) . (4.23) 
Due to Lemma fl. 21 for polynomials in z we have 

where Q is the matrix inverse to Sap and it is given by ()4.2|1 . Let r : — ^ Vm 
be the operator of contraction of a symmetric tensor with the matrix Q: 

By using formula ()4.23p and performing straightforward computations, we get 

A, = Q - ai - 2 r (C(2)) + r -\t' (%)) + ^ (S(3) V 5(3)) 
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A, 



Q 



a2-T (ai(2)) + 2 aiT (C(2)) - air (77(2)) 



1 1 
+ - aiT^ (%) - (C(4)) + - (7/(4)) + (^(3) V ai(i)) 

+ 6 (C(2) V C(2)) + ^ (r/(2) V r/(2)) - 3 (((2) V r/(2)) 

- ^ (%)) + ^ (%) V C(2)) -Ir' (5(4) V r/(2)) 

5 10 5 

- (^(3) V 5(3)) + y (5(3) V C(3)) - 3 (5(3) V r/o)) 

7 3^ 
+ Y^^' (%)V5(3)) + -r^ (%V5(4)) 

- y (-^O) V 5(3) V C(2)) + (% V 5(3) V r/(2)) 

(5(4) V 5(3) V 5(3)) 
+ 1^^' (5(3) V 5(3) V 5(3) V 5(3)) 



All terms in these expressions were evaluated above. Their substitu- 
tion gives expressions for the computed coefficients in terms of the curvature, 
the differential (i$ of the mapping $, and their covariant derivatives. 

Lemma 4.2. Let conditions of Theorem \4.1\ be satisfied. Then 



H. 



(2) 



Q'^^Q^'Q^'^Q^' + hI 



(2) 

(al3'y5e(r)eM:\) 
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where 



111 1 

72 ^ + 30 180 + 180 ' 

'^^ ^ = RRnR + -En-/? + i ~ R un — — 

36 ^ ^ V6 1^ 



3 



11 ^ 

17(2) ^ p p -'- p p/i 1/ p 

-"a/375 ~ g ^apn^tS J^fiaulS-n 7 <5 2Q 

- R^aRu/S + - RRfMau/S — Rfia;u/3 + 3 -F)ti/a/3^ 
+ 1— i? Be + — B B^J + — i? t \1/'^«\1>^ \1>'';; 

I 2^ -'•'/xiz-'i'gTr 2Q J-^efiTi/J^i' 20 ^^''^'^ 1 " ^ 

+ ^ -R;/x^'^a^/375 + ^2 -^M^^;" ~ 2 '^^"''^ ~ 2 '^'^"''^^ ^'^/3^'^75 
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5 2 1 \ 

Y2 Rpi-oiupRs^-K ~l~ 2 Rg/iauR 5/37r ~l~ g RfJ.auP;^Tr 1 ^ 7^ 5^ 



m(-R --R I M/^flVl/'' \l/;r ^ 

/5 15 \ 

+ ( ^ R^^^a^■,0 + ^ i??aM/3;^ ) ^''7^''<5^^< 

3 3 \ 

3 5 5 

4 iz o 

/5 3 

/in ^ \ ^ 



o 4 
3 \ 1 1 



5 5 

35 35 
35 35 

35 35 

- — R^,aup■,i^%'^''s^K'^nC^%e - —Rf^a'^'^p^-yse^cve 

/35 35 
/35 

+ i^M^^'/S ( y R^aup'^^^'^^'^s'^eCve + 35 ^^^^5e^Cr?e 

35 35 

- — Rc,f3-^^^-^^'se'^u('^%e + -TT Rf^a'^^is'^^^-^'^ss'^cve 
o Z 

35 / 

+ y + y *^*^^%*«cr?e 

35 35 

3^ 3^ 

35 35 

+ — R^^^\^''p^^^^^5e'^7r('^''rje + ^ ^a/37<5^<r,e 
O O 

35 35 35 

+ y ^'a/375e^Cr?e - y V^C^y^ " y ^';.a^'''/37fe^'Cr?^' 

35 35 35 

- -r ^a/375e*<*"^e - — ^'^a/3*''75*eC'?e - — ^'Ma^"'/375*eCr?e 
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35 35 



105 

35 



35 



^(2) 



315 



3 

3 

+ 2 ^'^^^^'^^^^'fe^^'^^^'^^e^A + ^^^a^'''/37^'^.5^''^<^'?e,.A 

3 



H. 



(2) 



3465 



+ 6 ^p,„^''/37^^5^%C^r,0,<^AV... 
- 4 ^f^„^"^^^^'^5^%^^^^^^e,.^A^^ 



g = |det (/ - d^{xo))\~^ 



and Ea/3, F^p^s are given by ()4.2ip . ()4.22p respectively. 



4.3 Isometric Case 

If the mapping $ is an isometry, all derivatives of (i$ vanish (see, 
example, ^3]) and from Lemma f4.2l we get the following lemma. 
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Lemma 4.3. Let conditions of Theorem \4.1\ be satisfied. Then 



Q 



'^^ 72 30 ''^ 180 ^ 180 



+ 



- ^ RR^p - ^ RR^a.'^^p + ^ i^i?^,^'^^^"^ 



a/3 



3 1 

o D 

1 1 



IMuaf}^ 7^ (5 



g(a/3g75) 



+ 
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g =\det[l - (i$(xo)) , 
and EajB, F^js^s o,re given by ()4.2ip . fl4.22p respectively. 

4.4 Comparison of the Results 

In this section we compare our expression for the coefficient Ai in a 
particular case when the mapping $ is an isometry 




(4.24) 



with the result obtained by Donnelly in jS]: 



= |det5| 



(4.25) 



where 



B = {I- d<i>) 



and all tensors are evaluated at the point Xq. 



From (jl^ it follows that 



ga/3 ^ ^a^^^/^M _ 



Since d^{xo) = and (/ — d^)B = I, we have 
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and (|4.24p takes the form 

= |det B| 1^ i? - ^ Ra(3B%B^^' - ^ - (5^)5°^ 

+ Ro.,puB^\B^\^\ {B\ - 5\) {B^^- - (J'^)-) 

After performing straightforward computations, we get the same expression 
as dOH). 



CHAPTER 5 



Conclusions 

Let us summarize the results obtained in this thesis. Our goal was 
to study the geometry of the fixed point set S of a smooth mapping $ on a 
Riemannian manifold (M, g) by computing the asymptotic expansion as t — * 0"^ 
of the trace of the deformed heat kernel U(t; x, of the Laplace operator 

on M. The fixed point set S was assumed to be a smooth compact submanifold 
of M. 

Multi-dimensional Gaussian integrals are expressed in terms of full 
contractions of the symmetrized products of symmetric tensors with a sym- 
metric (2,0)-tensor. In Chapter |21 we developed a technique for obtaining 
symmetrization-free formulas for such contracted products. The described al- 
gorithm can be easily implemented for automatic computations. 

In Chapter 01 we considered the case of a fiat manifold, M = M^. We 
showed the existence of the asymptotic expansion in the form 

/oo „ 
dvol(x)t/(t;x,<l>(x)) dvol(i/)afc(y), 

where m = dimS and are scalar invariants depending polynomially on 

derivatives of the function S{x) = cr(x, <l>(x)), and the inverse matrix of its 
Hessian matrix. We developed a generalized Laplace method for computing 
the coefficients and computed explicitly the coefficients Qq, ai, and a2 for 
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both zero- and one-dimensional fixed point sets S. We used the results from 
Chapter |^ to obtain a symmetrization-free form in the zero-dimensional case. 

Finally, in Chapter|3]we obtained the asymptotic expansion as t — 0"^ 
for an arbitrary curved manifold M and a zero- dimensional fixed point set 
S = {xq} in the form 

/oo 
dvol(a;) U{t; x, ~ ^ t^Ak , 

where Ak are scalar invariants depending polynomially on covariant deriva- 
tives of the curvature of the metric symmetrized covariant derivatives of the 
differential ci$ of the mapping $, and the matrix 

Q = {{I-d^fg{I-d^)]-\ 

evaluated at the point xq. The coefficients are expressed in terms of sym- 
metrized covariant derivatives of functions of the form /(x, We developed 
an algorithm for computation of such symmetrized covariant derivatives and 
realized it as the Python script listed in Appendix lA.il We computed explicitly 
the coefficients Aq, Ai, and A2. In a particular case, when the mapping $ is 
an isometry, our result for the coefficient Ai coincides with the former result 
of Donnelly [H]. 



APPENDIX A 



Computation of Symmetrized Covariant Derivatives 
A.l Generating Python Script: Source Code 

##################################################################### 

# Computes symmetrized covariant derivatives of the function 

# S(x)=\sigma(x,\Plii(x)) 

# on the fixed point set of the mapping \Phi(x). 
# 

# Input : 

# sSeparator - a separator which will be inserted between all terms 

# in the output 

# nDerivativeOrder - the order of the derivative to be computed 
# 

# Uses indices from IsAlphabet for derivatives and from 

# IsPrimeAlphabet for dummy indices in contractions 
# 

# Creates in the current directory file named "output.tex" which 

# contains symmetrized derivatives of the function S(x) up to 

# nDerivativeOrder written in LaTeX notation. 
# 

# Properties of the world functions, taken into account: 

# [\sigma]=0 

# [\sigma_\alpha] =0 

# [\sigma_{\alpha\beta\gamma}] =0 

# [\sigma_{ (\alpha_l\dots\alpha_n) }] =0 , n>=3 

# [\sigma_{\mu(\alpha_l\dots\alpha_n) ]■] =0 , n>=2 
# 

# An existing file with the same name will be overwritten. 
# 

# Copyright (c) Audrey Novoseltsev, 2005 

##################################################################### 

import OS 

from array import array 
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lsAlphabet=[ 
r"\alpha" , 
r"\beta" , 
r"\gamma" , 
r"\delta" , 
r"\eps" , 
r"\zeta" , 
r"\eta", 
r"\theta" , 
r"\kappa" , 
r"\lambda"] 

lsPrimeAlphabet= [ 

r"\mu" , 

r"\nu" , 

r"\xi", 

r"\pi" , 

r"\rho", 

r"\tau", 

r"\phi", 

r"\chi" , 

r"\psi", 

r"\oniega"] 

sSigma=r"\sigma" 
sPsi=r"\Psi" 

sSeparator="\n"+r"\displaybreak[0] \\ &"+"\n" 

iiDerivative0rder=6 

def NextDistribution(n,m,k,aI) : 
if al[0]==n-m-k+l: 

return False 
aIOold=aI [0] 
al [0]=1 

for i in range (l,k): 
if al[i] !=1: 

al[i]=al[i]-l 
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aI[i-l]=aIOold+l 
return True 



# Computing terms coefficients 
dCoef f icients={} 
dCoeff icients [0,0, ()]=! 

for n in range (1 ,nDerivativeOrder+l) : 
for m in range (n): 

for k in range (1 ,n-m+l) : 
al=array ( 'H' ) 
for i in ranged, k): 

al . append (1) 
al . append (n-m-k+1) 

while 1==1: 

#Compute the current coefficient 
coef =0 
if m!=0: 

coef =coef+dCoeff icients [m-1 ,k,tuple(al)] 
if al[k-l]==l: 

coef =coef+dCoeff icients [m,k-l , tuple (al) [0 :k-l] ] 
for i in range(k): 

if al[i]>l: 

al[i]=al[i]-l 

coef =coef+dCoeff icients [m,k, tuple (al)] 

al[i]=al[i]+l 
dCoeff icients [m,k, tuple (al)] =coef 
#Proceed to the next coefficient 
if not NextDistribution(n,m,k, al) : break 
dCoeff icients [n,0, ()]=! 



# Constructing LaTeX expressions 
f Output=f ile ( ' output .tex' , ' w+' ) 

f Output. wr it e( 
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"S=0" 

A-^^\r\°/°i°i°i°i "/ "/ "/ v v v v v v v v v v V v V V y y y y v v v y y y y v v v v y y y y y y y \ n " 

\Il/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/^ 

+"S_{"+lsAlphabet[0]+"}=0" 

-1- " \ T-,"/ y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y y \ n " 

\Il/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/^ 

+"S_{("+lsAlphabet[0]+lsAlphabet[l] +")}=" 
+"g_{"+lsAlphabet [0] +lsAlphabet [1] +"}" 
+"-2\ , "+sPsi+"_{ ( "+lsAlphabet [0] +lsAlphabet [1] +" ) }" 
+"+"+sPsi+"_{"+lsPrimeAlphabet [0] +" ("+lsAlphabet [0] +"} " 
+sPsi+""{"+lsPrinieAlphabet [0] +"}0_{"+lsAlphabet [1] +") }") 

for n in range (3 ,nDerivativeOrder+l) : 

f Output . write ( " \n%%%%%%%%%yoyoyo%%%%y,y,y,y,%%%%%%n%%%y,y,y,%%%%yoyoyoyo%%\n" ) 

f Output .write ("S_{(") 

for i in range (n): f Output .write (IsAlphabet [i] ) 
f Output . write ( " ) >=\n" ) 

nNumberOf Terms=0 

for t in range (2 ,n+l) : 

for m in range (t-1 , -1 , -1) : 
k=t-m 

al=array ( 'H' ) 

for i in ranged, k): 

al . append (1) 
al . append (n-m-k+ 1 ) 

while 1==1: 

#Check if the current term should be written 
if t==3 or (k==l and m>=2) or (k>=n-l and k>=3) : 
if not NextDistribution(n,m,k,aI) : 

break 
else : 

continue 
#Write the current term 

if nNumberOf Terms ! =0 : f Output .write (sSeparator) 
if k==l and m==l : 

fOutput . writeC'-") 
elif nNumberOf Terms ! =0 : 

f Output .write ("+") 
nNumberOf Terms=nNumber Of Terms+1 
if dCoeff icients [m,k,tuple(al)] !=1: 
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f Output .write(str(dCoeff icients [ni,k,tuple(al)] ) 
+"\, ") 
if t==2: 

fOutput .write ("g_{") 
for i in range (k): 

f Output . write (IsPrimeAlphabet [i] ) 
if m!=0: f Output .write (" (") 

for i in range (m): f Output .write (IsAlphabet [i] ) 
f Output .write ("} ") 
else : 

f Output .write (" ["+sSigma+"_{") 
for i in range(k): 

f Output . write (IsPrimeAlphabet [i] +" ' " ) 
if ni!=0: f Output .write (" (") 

for i in range(m): f Output .write (IsAlphabet [i] ) 

f Output .write ("}] ") 
f Output .write (sPsi+"~{"+lsPrimeAlphabet [0] +">{>_{") 
if m==0: f Output .writeC (") 
for i in range (m,m+al [0] ) : 

f Output . write (IsAlphabet [i] ) 
nextlndex=m+al [0] 
for i in range (l,k): 

f Output .write ("} "+sPsi+"~{"+lsPrimeAlphabet [i] 
+">{>_{") 

for j in range (next Index, next Index+al [i] ) : 

f Output .write (IsAlphabet [j] ) 
nextlndex=nextlndex+al [i] 
f Output .writeC')}") 
#Proceed to the next coefficient 
if not NextDistribution(n,m,k, al) : break 
f Output .write ("\n°/„ Number of Terms: "+str(nNumberOf Terms) ) 

f Output . close () 

A. 2 Generating Python Script: Output 

S = 0, 

Sa = 0, 
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7<5e) 

+ 10 [v.'(a/3]*''7^"5^) 
+ 20 [fr^V'(a/3]^''75^''s) 

+ 10[a^v'e(a]^^*%*^.) 

+ 10 [v.'(a/37]^''5^%) 

+ 10 [Vv'«'(a/3]*%*'5*^.) > 
'S'(Q/375eC) = — 6 5'^(Q,\[^^^^5e^) 

+ 10 ^^^^'^(^^^^'^5,^) 

+ 15 [ff;.V'(a/3]^''7^''5<) 
+ 45 [a^V'(a/3]*%5^%C) 

+ 45 [a^v'(a/3]^V^''c) 

+ 24 [v^'e(a]^''/37^''5^^eC) 



90 



+ 48 [v.'C'(a]^V^'^5.^^c) 
+ 36 [v.'«'(a]^V5^'.^^C) 

+ 3 [a^v'5'.']*^a*"/3*^75*%C) 

+ 3Kv'5'.']^''(a^"/3^S5e^"c) 
+ 8 Kv'5V']^''(a^"/37^^5e^"c) 

+ 10 [a^v'c%']*^a/3^%^^&^"c) 

+ 6[(T^V'5'.']*'^(a*"/375^l^*"c) 

+ 15 
+ 10 

+ 20 [a^V'(a/37]*''5*%C) 

+ 40 [a^V'(a/37]^^e*''c) 

+ 15 [a^V'e(a/3]^'^7*''5*^eC) 

+ 30 [v.'e(a/3]^%^"fc^^C) 

+ 45 [a^V'?'(a/3]*^7<5*''s^^C) 

+ 6[(T^V'5V'(a]^^*%*^5*%C) 

+ 12[(T^,,,5,,,(J*^^*%*«fe*'^C) 

+ 18 [a^...e-'(c.]^^^%*^.^"c) 

+ 24 [a^v'c'.'(a]^V^'5^^.^"c) 
+ 15[a^v.(a/375]*^*''c) 
+ 20 [v.'c'M7]*^^''.^^C) 
+ 15 [v.'5'.'M*%*"5^«.*"c) • 
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A. 3 Coincidence Limits of the World Function 
A. 3.1 Forth order 

111 

1 1 

g i^Rau^^ ~l~ Ra^uS^ ~\~ Rau^^) ~l~ 2 RcifJiS,!^ 

1 2 

~ 3 ~^ Ra^i(^u) 2 Ra{^i\i\y) 1 

2 

2 



A.3.2 Fifth order 

2 1 / 

+ 2 [(^a^l'u'■K'^,'] — - [{o'aa'R"' /xVvr');?'] 
1 1 

2 12 
--iff -Iff -Iff -2p 

= g RMa\u)[C,n) - g %|a|7r)(M;i^) , 
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2 5 5 

RaiP{i,;u) - Ra/3afi;u) 

Y %ICIi.){q;/3) - — R{^,\c.\u)^3■,^, - — R{a\m(p;^) ■ 



A.3.3 Sixth order 

[(^^l'u'{af3'y5)\ = — [cr^j'(a/37<5)i/] , 

[0"/x'a(/375i/)] = Y ([c"At'Q{/37(5)!^] + [o'fj'a(/37|!y|(5)] + n> a{p\u\-i&)\ + [o"fj'Qi^(/375)] ) 

= Y [V"{/37<5)i'] + Y [V"(/37k|5)] + 2 [V"{/3kl75)] 

- Y [{(y^l'aR''a{P\u\)■a5)\ 

13 1 

= Y K'q(/375)i.] + Y [o-^'q(/37|i.|5)] " 2 [(o"M'«(/3^"|a|7kl);5)] 

la 1 a 1 

- - [{(y^l'aaR''{f^-1\u\)■,5)\ - Y [0-;,'a{/37]-R"|a|5)i/ + -RiJia{l3\u\-a&) 

3 3 

= [<7Ai'a(/375)i/] — Y [0"A''a(/37]-^"|a|'5)i^ ~ ^ ^^' aa{fi\R"' 

3 1 

- Y [(y^L'a{l3\a\]R''-,5)v - - [(^ ^l' a(p-i]R"' \a\5)v 

- 2 V'' ^l' aa{p]R'^ -y&^u ~ ^ [c"^j'a(/37] -R*^ |a|(5)v + ^ -R/xa(/3|y| ;7(5) 

3 5 

= [(^ ^^' a{f3^5)u\ — 2 ['^M'a(/37]-^"|a|'5)!^ ~ ^ [o"/i'aa(/3]-R"75)i/ 

3 1 

- Y [Va(/3k|]-^''7<5)f + Y R^lOL{|3\ul-|5) , 

3 5 

[0"/x'«(/375)v] = [cr^'a(/375i/)] + 2 [^/i'a(/37]-R'^|a|5)i^ + ^ [o"M'aa{/3]-^"75)'^ 

3 1 

+ Y [V"(/3|a|]-R"75);^ - Y Ri,a{(i\u\-n&) 
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3 

- R^li(3\a\rM) - 


7 


- 2 -RM/3|a|7^''|a|5), 


p pa 


5 


p pa 

O 


~ g -RM{/3|aa|-R"7(5)!^ " 


^ P pa 


g -R^to{/3|a|-R''75)!/ 


- 1 -RMa(/3|i.|;75) 






3p 




-R/x(/3|a|7-R°5|a|j/) 


- - -RM/3|a|7-^"l"|5)!^ 


p pa 

2 ^^aa{l3^ 7(5)i/ 


" 2 RfJ.(/3\aa\R"''yS)u 


p pa 


3 

o Rlia{l3\a\R ■y5)u , 





3 17 

[0"/,'(q/375)z.] = -— ^/.{Q/3|i.|;75) " ^ -R/.(a/3|i.|;75) " — -RM(a|a|/3-R%5)v 

115 

J? I?a D Da p Da 

— 2 -"-/^{"lal/?-"- 7(5)!^ ~ 2 'tt^J.(a|3\a\'t^ ■yS)u — g -rt^(Q|a|/3-rt 7.5)1/ 

3 

~ o Rfi{af3\a\R ■y5)iy 
O 

2 149 17 

= -RM(a/3k|;7<5) - Y;^Riiia\a\pR°''yS)u " ^ -R^Q/Slal 

2 8 

= ^ ^M{a|i'|/3;75) + Raia^pR" ^\u\5) , 
0'AtV?'(a/37)] = [crMV(a/37)];5 — W ^' u' (a/S-y)^] 



[V(«/375)i'] - 4 ([V(a/37)5'^] + ['^/^'(a/3|5|7)i'] + [V("ISI/37)f ] + [ V5("/37)i'] ) 

111 

= ^ [0"M'(a/37)?'^] + ^ ['^M'(a/3|Cl7)'^] + g [^^''("1^1/57)'^] 

13 1 

= ^ [0"m'(«/37)C'^] + ^ ['^M'(a/3|Cl7)'^] " 2 [(^M'a-R"(a/3|€|);7)i'] 

3 

= [crM'(o/37)ef] - 2 

- 2 ['^M'a(a|i/|]-R''/37)|C + 2 -^MQ/3|el;7)i^ 

= [0"^'(a/37)5i/] - 2 [0-^'a(a|,.|]-R''/37)g + - -R^{a/3|e|;7)i/ , 
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2 8 

= ^ Rfiia\u\l3;'ri) + ^a(a|/.|/3^''7|^|5) 

4 ^1 

5 11 

5 2 8 

- — R{aMm;'r)u + - Rfi{a\u\P;'ri) + ^a(a|/.|/3^''7|i.||) 

4 a 1 

+ 2 ^(M|a|;.)(a^''/37)5 " " Rf,(am;i> 

= "2 ^M(c^kl/3;7)C + 5 ^m(«I'^I/3;7«) 

8 a 4 

+ — RaiaMf^R" -yliylO + g ^(M|a|;.)(a^''/37)e 

= "2 ^M(a|^l/3;7)C + Yo ^m(«I'^I/9;7)« + Yo ^M(a|^l/3;l5l7) 

+ ^M(c^k5l;/37) + Yo ^m5i'(c^;/37) + Y5 

2 2 2 

+ Y^ Ra{a\^l\|3R''\^u\^) + Y^ i?a(Q|MC^''/3k|7) + Y^ Ra^fJ-iaR'' p\u\y) 
4 

+ 2 ^(M|a|i.)(a-R''/37)? • 

Symmetrization in (/i, z/, ^) is assumed in all terms of the following expression, 

2 14 

2 2 2 

= i?M(c^kl/3;l?l7) + 5 ^a(a|i.|/3i?''|/.C|7) + ^ ^m«i^(«^''/3|5|7) 

1 4 ^ 

+ Y^ ^/^(a|H/3;l«l7) - ■^Ra^l{a\u\R''m\■y) 

3 4 ^ 

= "YO ^M"kl/3;ICl7) ~ 3 ^a/^lali/l^^/Sl^h) • 
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